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ABSTRACT 


The  similarity  laws  for  the  experimental  study  of 
aerothermoelastic  phenomena  associated  with  flight  at  hyper¬ 
sonic  Mach  numbers  are  derived.  This  study  considers  in 
detail  the  parameters  associated  with  the  external  pressure 
and  heat  flux  imposed  on  the  structure,  the  heat  transfer 
distribution  within  the  structure  (for  both  ablating  and 
non-ablating  materials)  and  the  stress-deflection  behavior 
under  short-time  and  long-time  loading  conditions.  The 
structural  parameters  associated  with  nonlinear  material 
properties  at  high  temperatures,  plasticity,  creep  and 
fatigue  are  briefly  considered. 

Similitude  for  the  general  aerothermoelastic 
problem  is  possible  only  for  a  scale  ratio  of  one  and  for  a 
thermal  and  pressure  environment  identical  with  the  prototype. 
The  basic  conflict  arises  when  mutually  satisfying  the  require¬ 
ments  of  same  temperature,  Reynolds  number,  surface  radia¬ 
tion  parameter  £oTo3l/k0  and  the  aeroelastic  parameter  U /eo  • 

Means  of  dealing  with  this  conflict  are  suggested. 
These  include  looking  at  more  specialized  situations  pecul¬ 
iar  to  both  blunt  bodies  and  affinely  related  slender  bodies 
in  continuum  and  free  molecule  flow,  the  use  of  "incomplete 
aerothermoelastic  testing"  in  which  pressure  and  thermal 
loadings  are  estimated  in  advance  and  applied  artificially 
to  the  model,  and  the  use  of  "restricted  purpose"  tests 
investigating  a  particular  facet  of  the  complete  aerothermo¬ 
elastic  problem. 

Finally,  a  survey  of  some  available  types  of  test 
facilities  appropriate  for  these  investigations  is  given. 
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1. 


INTRODUCTION 


The  present  report  constitutes  the  second  phase 
of  a  study  concerned  with  the  similarity  laws  required  for 
experimental  aerothermoelastic  studies.  The  first  phase 
was  primarily  concerned  with  flight  Mach  numbers  extending 
into  the  higher  supersonic  range  (M  <_  3.5)  and  with  maximum 
structural  temperatures  of  approximately  1000°F .  Materials 
with  elastic  stress- strain  properties  only  were  considered 
and  structural  radiation  was  assumed  negligible.  The 
area  of  interest  of  this  first  phase  was  restricted  to 
region  I  of  Figure  1  and  the  results  are  presented  in 
Ref.  1.  It  was  found  that  the  basic  similarity  parameters 
for  combined  aerodynamic,  heat  transfer  and  stress  deflec¬ 
tion  studies  are  the  free  stream  Reynolds  number  Re  m  > 

Mach  number  and  the  basic  aeroelastic  parameter  1*/E0  • 

In  the  present  report  the  similarity  requirements 
associated  with  Mach  numbers  extending  into  the  hypersonic 
range,  regions  II  and  III  of  Fig.  1,  are  considered. 
Essentially  the  basic  ideas  of  similarity  and  non- dimens ion - 
alizat ion  of  Ref.  1  are  utilized  and  extended  to  cover 
these  cases.  Some  of  the  assumptions  made  there  are  no 
longer  applicable.  The  higher  temperatures  associated 
with  hypersonic  flight  introduce  new  and  considerably  more 
complex  phenomena  relating  to  aerodynamic,  material,  and 
structural  characteristics.  The  similarity  laws  pertaining 
specifically  to  hypersonic  pressure  and  heat  transfer  simil¬ 
arity  must  now  be  investigated.  Such  additional  phenomena 
as  real  gas  effects,  structural  radiation,  material  ablation, 
nonlinear  behavior  of  material  properties  at  high  temper¬ 
atures,  plasticity,  creep  and  fatigue  are  also  to  be  con¬ 
sidered  now. 


1 


This  current  study  will  first  discuss  in  detail 
the  parameters  associated  with  the  aerodynamic  pressure 
distribution,  aerodynamic  heat  transfer,  heat  conduction 
for  non-ablating  and  ablating  bodies,  and  the  stress- 
deflection  behavior  at  high  speeds  and  temperatures.  The 
possibility  of  mutually  simulating  the  combined  aerodynamic, 
heat  transfer,  and  structural  behavior  problem  will  then  be 
discussed.  Again  as  in  Ref.  1,  some  basic  conflicts  will 
arise,  and  wherever  possible,  some  means  of  alleviating 
them  will  be  discussed.  Applications  will  then  be  made  to 
both  boost  glide  and  ballistic  type  vehicles.  Finally, 
some  of  the  current  research  facilities  available  for  testing 
these  aerothermoelastic  models  will  be  described. 
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2. 


AERODYNAMIC  SIMILARITY 


2.1  General  Features 

The  two  aerodynamic  quantities  of  primary 
importance  in  aerothermoelastic  modeling  are  the  pressure 
distribution  and  the  heat  transfer  rate  imposed  on  the 
structure  by  the  external  flow.  The  general  similarity 
parameters  relating  to  these  aerodynamic  quantities  were 
obtained  in  a  previous  report.  Ref.  1,  for  flight  extending 
into  the  supersonic  range  (M <  3.5).  It  was  shown  there 
that  if  the  following  non-dimensional  quantities 


f 


) 


cc  = 


A 


f^-ao 


(2.1) 


•7ere  introduced  into  the  general  equations  of  motion  of  a 
compressible,  viscous,  heat-conducting  perfect  gas,  the 
following  main  similarity  parameters  would  result. 


77 

Tgo 


same  geometry  ,  p  ,  /Cr ,  A 


(2.2) 


If  the  above  similarity  parameters  were  made  the  same  for 
both  model  and  prototype,  the  non-dimens ionalized  equations 
of  motion  would  be  identical  for  both,  and  would  therefore 
have  identical  solutions  for  the  unknowns  f,  j? >  7J  u,  , 
etc.  Hence  measurements  of  f  >  -p }  T >  ^  etc.  on  a 
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on  the 


model  can  then  be  used  to  find 
prototype  by  application  of  Eqs.  (2.1).  Therefore 
aerodynamic  similarity  in  both  pressure  distribution  and 
heat  transfer  will  prevail  if  the  similarity  parameters 
of  Eq.  (2.2)  are  maintained  equal  for  both  model  and 
prototype.  Further,  it  was  indicated  in  Ref.  1  that  M  , 
and  to  a  lesser  extent  y  ,  are  (provided  Re^  is  suffi¬ 
ciently  large)  the  principal  parameters  for  pressure  dis¬ 
tribution  similarity,  while  Re^  and  to  a  lesser  extent 
,  are  the  principal  heat  transfer  parameters. 

The  unsteady  term  UL°/l  and  the  reference  temperature 
ratio  11/ X»  >  serve  to  define  the  reference  time  t* 
and  the  reference  temperature  X  ,  respectively.  The 
remaining  parameters  j*  >  &  require  that  the 

variations  of  ,  /C -f}  A.  with  temperature  be  of  the 
same  form  for  model  and  prototype,  and  are  readily  met 
in  the  range  up  to  1000°F. 

At  much  higher  Mach  numbers,  when  one  gets  into 
the  hypersonic  speed  range,  an  important  new  property 
comes  into  play,  namely,  the  Mach  Number  Independence 
Principle  discussed  by  Hayes  and  Probstein  in  Ref.  2. 

This  states  that  at  a  sufficiently  large  free  stream 
Mach  number  the  flow  over  a  body  becomes  inde¬ 

pendent  of  M^  .  Considering  the  same  general 
equations  of  motion  of  a  compressible,  viscous,  heat- 
conducting,  perfect  gas  as  given  for  the  supersonic  case 
in  Ref.  1,  one  now  finds  it  convenient  to  match  conditions 
at  the  stagnation  point  and  to  non-dimensionalize  with 
respect  to  stagnation  point  rather  than  free  stream 
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quantities  (except  for  f>  and  u.  )  as  follows. 


p 

ft*  J 

ft 

,  T~ 

T 

To 

—  66 

66  =  -H-  j 

'  U 

i  > 

7-P 

) 

t 

fo 

(2.3) 
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A  • 

*rs  1 

A 
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TTs 

Introducing  these  into  the  equations  of  motion  as  was  done 
in  Ref.  1,  and  noting  the  following  stagnation  point 
approximations  for  high  /M<»  flow. 


'fs  ~ 


(r+£L 

z  (y*-0 
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p.  ~  p  (2.4) 

(t"i)(r+-')  /<2° 

J.S  =  T/t 

-r  ~  J/— 

S  2  f 

one  obtains  then  the  following  main  similarity  parameters 
for  the  hypersonic  flow  of  a  perfect  gas, 

Res  }  frs  >  Y  )  ,  same  geometry 

=  (2.5) 

A 

where  ^LJPs  and  /*.  =  /Cf’sY9/k$  • 


*The  factor  P  =  fo‘cpA (ts)  accounts  for  variations  in  .  For 
Xp  =  constant,  f,  =  l.  If  £p  is  simulated,  then  f,  ,  will 
also  be  the  same  for  both  model  and  prototype  in  the  above 
perfect  gas  equations. 
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Comparing  the  above  with  the  previous  similarity  parameters 
of  Eq.  (2.2),  it  is  seen  that  the  dependence  on  the  Mach 
number  Afc  has  been  dropped.  This  was  made  possible  by 
the  high  M &  approximations  of  Eq.  (2.4),  which  carry  the 
implication  that  the  free  stream  quantities  /fa  and  7^,  are 
very  much  smaller  than  the  corresponding  stagnation  point 
quantities  and  Ts  ,  and  hence  can  be  neglected  in 
comparison.  Also  it  is  seen  that  the  /t?  condition  for 
the  reference  temperature  %  replaces  the  previous  ~T° /T*> 
condition  of  Eq.  (2.2).  In  the  similarity  parameters  of 
Eq.  (2.5),  the  Rcs  again  is  associated  with  viscous  and 
heating  effects,  and  can  be  neglected  when  considering 
pressure  distribution  similarity  around  blunt  bodies.  For 
slender  bodies  however,  there  appears  to  be  a  viscous 
interaction  between  the  thickness  of  the  boundary  layer 
and  the  external  non-viscous  flow  at  these  high  Mach 
numbers,  and  the  Res  condition  does  play  a  role  in  the 
aerodynamic  pressure  distribution  similarity  here.  Of  the 
remaining  parameters,  it  is  to  be  noted  that  although  the 
Prandtl  number  fy-  can  be  assumed  fairly  constant  over  a 
wide  temperature  range,  the  simulation  of  <cf  ,  pi ,  Jf 
is  somewhat  more  difficult  over  these  extended  temperature 
ranges  associated  with  high  Mach  numbers.  As  mentioned  in 
Ref.  1  though,  any  gas  whose  yC-pj  and  properties  can 
be  represented  by  the  form  T“  will  satisfy  the  conditions 
on  jC. f ,  and  A.  .  Of  course  operating  with  the  same 

gases  at  the  same  temperatures  for  both  model  and  prototype 
will  also  insure  simulation  of  these  quantities. 

The  special  case  of  slender  affinely  related 
bodies  will  be  treated  more  in  detail  subsequently  in 
Section  2.2.  It  will  turn  out  there  that  further  relax¬ 
ations  in  similarity  can  be  made  for  these  slender  bodies 
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provided  the  parameter  7  is  large*  where  T  is  the 
thickness  ratio  of  the  profile. 

Thus,  for  general  hypersonic  similitude  of  a 
perfect  gas,  a  relaxation  of  the  Mach  number  condition  is 
possible  if  is  sufficiently  large.  It  is  interesting 

to  notice  that  for  very  low  speeds  (incompressible  flow)  as 
well  as  for  these  hypersonic  speeds,  the  Mach  number  again 
ceases  to  be  an  important  similarity  parameter  in  aero¬ 
dynamic  flow. 

A  survey  of  hypersonic  flow  problems  is  given 
in  Refs.  3  and  4.  Aside  from  the  Mach  Number  Independence 
Principle,  the  distinguishing  features  of  this  flow  regime 
are  the  closeness  of  the  leading  edge  shock  wave  to  the 
body  surface  and  the  attendant  high  temperatures  which  may 
result  in  real  gas  phenomena.  The  closeness  of  the  shock 
wave  to  the  surface  has  enabled  the  flow  to  be  treated  as 
Newtonian,  particularly  over  a  blunt  nosed  body,  in  which 
case  the  flow  characteristics  become  independent  of  . 

This  is  another  manifestation  of  the  Mach  Number  Indepen¬ 
dence  Principle.  For  a  slender  body  this  closeness  of  the 
shock  wave  may  result  in  an  interaction  with  the  viscous 
boundary  layer  inducing  higher  pressures  and  heat  transfer 
rates,  and  requiring  consideration  of  Reynolds  number  Re$ 
effects  on  the  pressure  distribution. 

The  real  gas  phenomena  assume  importance  at  the 
high  temperatures  associated  with  >  7  (for  air),  and 
involve  consideration  of  dissociation  and  recombination 
processes  and  their  associated  relaxation  times.  They 
affect  mainly  the  heat  transfer  rather  than  the  pressure 
distribution  similarity.  To  simulate  these  real  gas  effects, 
it  is  essential  at  least  to  maintain  identical  temperatures 
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in  identical  gases  over  model  and  prototype.  From 

Eqs.  (2.3)  and  (2.4) ,  this  requires  the  free  stream 
velocity  V  to  be  maintained  the  same  for  both  model  and 
prototype,  and  hence  the  parameter  U  is  added  to  the 
group  given  in  Eq.  (2.5).  For  more  accurate  simulation  of 
real  gas  effects,  it  is  desirable  to  maintain  identical 
pressures  as  well  as  identical  temperatures  over  model  and 
prototype,  thus  assuring  that  the  gas  is  in  the  same  state 
for  model  and  prototype.  It  is  then  seen  from  Eqs.  (2.3) 
and  (2.4),  that  the  free  stream  density  as  well  as  U 
should  be  maintained  similar  for  model  and  prototype.  This, 
however,  causes  difficulty  in  fulfilling  the  Reynolds 
number  similarity  condition. 

The  similarity  parameters  associated  with  hyper¬ 
sonic  phenomena  will  be  considered  in  more  detail  in  the 
following  sections  with  regard  to  aerodynamic  pressure  and 
heat  transfer  similitude  individually.  This  separation 
is  made  possible  by  the  Prandtl  boundary  layer  concept.  By 
looking  at  specialized  situations,  it  will  be  seen  that  it 
is  possible  to  provide  further  relaxations  0f  the  general 
conditions  for  supersonic  flow,  Eq.  (2.2),  and  of  the 
general  conditions  for  hypersonic  flow,  Eq.  (2.5). 
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2.2  Aerodynamic  Pressure  Distribution  Similarity 
a)  Blunt  Bodies 


At  high  Mach  numbers,  the  flow  over  the  nose  of 
a  blunt  body  may  be  treated  as  Newtonian  and  the  pressure 
expressed  as  (Ref.  3) 

X  ~  sin* ft  (2.6) 

-Ts 

where  ft  is  the  angle  between  U  and  a  surface  element  (local 
angle  of  attack)  .  For  M#,  greater  than  about  3, 


(2.7) 


(r+o 


Tests  reported  in  the  literature  have  shown  that  this 
insensitivity  to  Mach  number  is  valid  for  as  low  as  3. 

Thus  the  pressure  distribution  over  the  nose  of  a  blunt 
body  may  be  adequately  simulated  by  reproducing  only  the 
same  geometry  ( ft  ) ,  provided  both  model  and  prototype 
operate  above  M  3  .  The  above  relations  also  follow 

from  application  of  the  Mach  Number  Independence  Principle 
which  led  to  the  parameters  given  in  Eq.  (2.5).  Since 
the  viscous  boundary  layer  does  not  affect  the  pressure 
distribution  in  the  neighborhood  of  such  blunt  bodies,  the 
Res  and  Fys  conditions  can  be  disregarded.  The  y  condi¬ 
tion  is  here  absorbed  in  the  definition  of  the  stagnation 
pressure,  while  the  l/i* /L  condition  can  be  neglected  if  the 
flow  is  assumed  steady.  Further,  the  pressure  distribution 


^L. 

ft  includes  the  instantaneous  velocity  of  the  boundary  as 
well  as  the  instantaneous  slope  in  its  definition. 
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is  not  much  affected  by  real  gas  phenomena  here--any  real 
gas  effects  are  manifested  in  a  decrease  in  the  value  of  y 
(from  1.4  to  about  1.2  for  air)  which  alters  the  stagnation 
pressure  Eq.  (2.7)  only  slightly. 

b)  Sharp-Nosed  Slender  Bodies 

The  hypersonic  similarity  rules  for  sharp  nosed 
slender  bodies  of  similar  shape  in  inviscid  flow  have  been 
developed  in  Refs.  5  and  6.  For  a  family  of  affinely 
related  slender  bodies  of  thickness  ratio  T  in  an 
inviscid  perfect  gas  hypersonic  flow.,  the  pressure  may  be 
expressed  as 

^  r"ir(y/L  ;  r,  y,  ,  uHl)  (2.8) 

where  Tf  is  a  function  of  the  non-dimensional  parameters 
shown.  In  the  above,  Uo  refers  to  any  deflection  of  the 
body  and  U°/lt  requires  this  to  be  in  the  same  ratio  as  the 
thickness.  The  unsteady  term  Ui*/L  above  defines  the 
reference  time  to  for  unsteady  flow  phenomena,  and  enters 
through  the  boundary  conditions  on  the  displacement.  Eq.  (2.8) 
can  be  rewritten  in  the  form. 


r(n-v) 


X  *.  TT  (*//.-,  A7.r,  r,&>  <£■)  <2-»: 


From  the  above  it  is  seen  that  the  non-dimensional  parameters 
required  for  similarity  of  hypersonic  flows  are. 


M„r,  tr ,  s  ,  ^  ,  uu 


(2.10) 
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where  8  represents  the  body  thickness  distribution  of 
the  family  of  bodies  with  different  r  ,  and  insures  that 
they  be  affine.  The  above  are  applicable  for  M L  as  low 
as  2.5. 

At  very  high  Mach  numbers  such  that  M&T  1 
even  though  TCI  ,  the  function  TT  becomes  independent 
of  /HbT  and  approaches  a  limit.  Also,  the  term 

becomes  small  in  comparison  with  TT  and  hence  can  be 
neglected.  The  N\«T  condition  then  drops  out  of  Eq.  (2.10) 
leaving  only 


r  U» 

Y  >  ^  >  LT  '  L 


(2.11) 


as  the  required  similarity  when  AA^T  >>  1  •  The  above 
is  a  more  relaxed  form  of  the  Mach  Number  Independence 
Principle,  Eq.  (2.5),  for  no  viscous  effects,  since  now 
the  thickness  ratio  T  need  not  be  simulated,  only  its 
distribution  8  .  The  presence  of  real  gas  effects  would 

result  mainly  in  a  somewhat  lower  value  of  y  here,  but 
would  probably  not  affect  the  inviscid  pressure  distribution 
marked ly . 

At  hypersonic  speeds,  the  shock  wave  about  these 
slender  bodies  comes  close  to  the  surface,  and  causes 
interaction  with  a  relatively  thick  viscous  boundary  layer. 
Significant  changes  in  the  pressure  distribution  may  occur 
as  a  result.  This  viscous  boundary  layer,  which  is  ordin¬ 
arily  neglected  in  determining  the  aerodynamic  pressure 
distribution,  must  now  be  considered  for  these  slender 
bodies.  It  brings  in  additional  Reynolds  number  effects 
into  the  simple  inviscid  similarities  of  Eqs.  (2.10)  and 
(2.11). 
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This  viscous  interaction  problem  is  considered 
by  Hayes  and  Probstein  in  Ref.  2.  The  flow  field  outside 
the  boundary  layer  on  a  slender  body  is  considered  to  be 
the  same  as  an  inviscid  flow  field  about  the  same  body 
increased  in  its  dimensions  by  the  boundary  layer  displace¬ 
ment  thickness.  It  is  important  then  that  the  total  thick¬ 
ness  distribution  including  the  displacement  thickness  be 
affinely  related.  This  requires  that  the  ratio  of  boundary 
layer  displacement  thickness  &  to  body  thickness  rL 
be  maintained  the  same,  thus  introducing  the  additional 
similarity  parameter,  ,  to  those  of  Eqs .  (2.10)  and 

(2.11).  An  analysis  of  the  hypersonic  boundary  layer  in 
Ref.  2  reveals  that  the  parameter  S  j^L  can  be  replaced  by 


r 


■  FTTul 

v 


(2.12) 


provided  that  is  the  same,  and  that  the  non-dimensional 

temperature  distribution  ~7w/~j~s  at  the  wall  surface  of  the 
body  be  maintained  the  same  on  both  model  and  prototype.  In 
the  above,  fw  and  /*w  refer  to  some  point  at  the  wall 
surface.  Alternatively,  the  parameter  X  can  be  combined 
with  the  hypersonic  similarity  parameter  M<x>  T  of  Eq.  (2.10) 

•k 

which  after  some  manipulation  gives  the  alternate  parameter. 


X 


(2.13) 


r  -  i  -  ^ 

"it  can  be  shown  that  c<  T  foo  and  that  lw  ct  Mco  Xo 

Also  sfC  oC  where  s  is  the  exponent  in  the  approximation 

«C  T  $ 
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where 


C  =  p*  X, 


Adding  these  terms  to  Eq.  (2.10),  the  pressure 
distribution  similarity  for  a  family  of  affine  bodies  in 
hypersonic  flow,  including  the  viscous  interaction  effects 
requires  the  following  parameters  for  low  and  moderate 

M„r ' a  , 


,  jr 


Uio  MjF  tZj 

L  '  7*r  ;  * 


(2.14) 


For  very  high  Mach  numbers  where  Moo'Y 
meters  become. 


the  para- 


Y ' 


Uo 

ur 


U-to  -i— _ 

~lT  '  r  [HYUl 

i 


(2.15) 


It  is  interesting  to  note  that  the  above  is  a  more  relaxed 
version  of  the  general  Mach  Number  Independence  Principle, 
Eq.  (2.5),  since  Res  now  appears  in  combination  with  T 

rather  than  separately. 

In  the  boundary  layer  theory  used  to  obtain  the 
similarity  parameters  of  Eq.  (2.14)  and  (2.15),  it  was 
assumed  that  the  boundary  layer  is  a  perfect  gas  of  constant 
y  }  but  that  the  viscosity may  be  variable.  To  account 

*The  fifth  parameter  of  Eq.  (2.15)  as  well  as  the  sixth  of 
Eq.  (2.14)  can  be  written  alternatively  as 
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more  accurately  for  real  gas  effects,  it  is  necessary  to 
use  the  same  gas  and  to  maintain  at  least  identical  temper¬ 
atures  within  the  boundary  layers  of  both  model  and  proto¬ 
type.  This  is  achieved  by  additionally  maintaining  7?  the 
same,  which  from  Eq.  (2.4)  implies  U  to  be  the  same  (for 
same  gas).  For  the  low  and  moderate  case,  the 

similarity  parameters  of  Eq.  (2.14)  then  increase  to 

$ ,  -jyy  i  -j—  t  is  ,  T  f^L ,  Ik_  ,  same  gas  (2.16) 

The  same  relations  hold  without  the  term  when  M*>7  >yl.. 

These  constitute  Similarity  Cases  (B)  of  Ref.  2. 

For  even  closer  representation  of  the  real  gas 
effects  in  the  boundary  layei;itis  also  proposed  in  Ref.  2 
to  maintain  identical  pressures  as  well  as  identical  temper¬ 
atures  over  both  model  and  prototype  thereby  assuring  that 
the  gas  is  in  the  same  state.  This  may  be  achieved  by 
maintaining  the  same  f>  u'1' as  well  as  the  same  Ts  (hence 
C 7  )  on  both  model  and  prototype  (see  Eq.  2.9) .  Intro¬ 

ducing  these  into  the  low  and  moderate  Y  case,  the 
similarity  parameters  of  Eq.  (2.16)  increase  to 


M.t.  S, 


HQ 

L'T 


Iz.  ,  V 


Z- 

'YL. 


' r5 


same  gas 


(2.17) 


The  same  relations  without  the  M*,?  term  hold  for  the 
Mao'*  >>  ±  case.  These  constitute  Similarity  Cases  (A)  of 
Ref.  2. 
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It  is  probably  sufficient  for  many  cases  of 
aerodynamic  pressure  distribution  similarity  to  utilize 
Eqs.  (2.14)  and  (2.15).  It  is  to  be  noted  that  for  small 
values  of  the  parameter  X  of  Eq.  (2.13),  viscous  inter¬ 
action  effects  are  weak  and  can  be  neglected  thus  dropping 
out  the  DC  and  associated  temperature  effects  and  leading 
back  to  the  inviscid  pressure  distribution  similarity  of 
Eqs.  (2.1o)  and  (2.11). 

c)  Blunt  Nosed  Slender  Bodies 

The  leading  edges  of  slender  bodies  are  generally 
blunted  to  reduce  the  high  heat  transfer  rates  at  supersonic 
and  hypersonic  speeds.  It  has  been  found  experimentally 
that  even  slight  blunting  may  induce  significant  inviscid 
pressure  distributions  over  a  body  surface  at  hypersonic 
speeds, (Ref.  7).  On  the  basis  of  hypersonic  small  pertur¬ 
bation  theory,  Cheng  in  Ref.  8  derives  the  following  addi¬ 
tional  similarity  parameter  to  account  for  these  blunt  nose 
effects  on  pressure  distributions. 


mJ  a  j- 


for  plane  surfaces 


for  bodies  of  revolution 


(2.18) 


In  the  above,  d  is  some  typical  dimension  of  the  blunt  nose 
(the  diameter  for  cylinders  and  spheres)  and  k  is  the  nose 
drag  coefficient,  which  varies  between  1,  4/3,  2  for  spheri¬ 
cal,  cylindrical,  and  flat  noses,  respectively,  in  hyper¬ 
sonic  flows.  Since  k  is  not  particularly  sensitive  to  the 
geometry  and  since  presumably  both  model  and  prototype  will 
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have  the  same  general  form  of  blunting,  the  dependence  on  k 
in  the  similarity  parameter  of  Eq.  (2.18)  can  be  dropped. 

The  similarity  parameter  of  Eq.  (2.18)  has 
correlated  data  for  5"  2.  2.  20  quite  reasonably  for 

the  -f/Tco  on  a  flat  plate  afterbody  for  different  combinations 
of  Woo  and  d,(Ref.  8).  The  similarity  parameter  of  Eq.  (2.18) 
also  occurs  in  the  "blast-wave  solution"  for  blunt  nosed 
bodies  (Ref.  9). 

The  basic  similarity  parameters  for  aerodynamic 
pressure  distribution  on  blunt-nosed  slender  bodies  in 
hypersonic  flow  can  be  obtained  by  adding  the  bluntness 
parameter,  Eq.  (2.18),  to  those  of  Eq.  (2.14).  For  lifting 
surfaces,  these  become. 


Af,r ,  r 


s  , 


o 

Lr 


Ubo 

L 


Ik. 

Ts 


K-  ,  Wj  A  J- 


(2.19) 


The  above  presumably  account  for  both  viscous  interaction 
and  blunting  effects.  Again,  for  >>  1  ,  the  depen¬ 

dence  on  Mo,'}'  can  be  neglected.  Also  the  above  can  be 
made  more  precise  by  using  the  conditions  of  Eqs.  (2.16) 
and  (2.17).  In  maintaining  similarity,  the  blunting  effects 
are  controlled  by  the  independent  quantity  d  which  can  be 
varied  at  one's  disposal.  It  is  to  be  noted  that  at  small 
values  of  the  parameter  A/}J  A  ^  ,  the  blunting  effects 

become  small  and  can  be  neglected. 

d)  Free  Molecule  Flow  Regime 

At  extremely  high  altitudes  the  air  may  be 
rarefied  to  the  extent  where  continuum  concepts  no  longer 
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apply.  The  parameter  distinguishing  flow  regions  is  the 
Knudsen  number  Kn  —  ^ /L  >  where  ^  is  the  molecular  mean 
free  path  and  L  a  characteristic  body  dimension.  In  terms 
of  free  stream  flow  paramters, 

oC  MfL.  (2.20) 

Free  molecule  flow  is  considered  to  be  established 
when  y  IO  >  according  to  the  criterion  established 

by  Tsien,  (Ref.  10)  Reference  11  indicates  that  free  mol¬ 
ecule  flow  in  the  case  of  a  cylinder  can  also  be  character¬ 
ized  by  Kh  >  2  .  Generally,  for  satellite  velocity 

vehicles,  the  free  molecule  regime  can  be  considered  as 
starting  above  about  90  miles  altitude. 

Utilizing  concepts  from  the  kinetic  theory  of 
gases,  the  pressure  distribution  about  a  body  in  free 
molecule  flow  is,  (Ref.  12) 


'/Z  /U  u/sl 


(z-f)  s  SVA'  0 

/tT 


(2.21) 


■h 


z 

■h  S 


where  f  is  the  fraction  of  molecules  that  are  diffusely 
reflected  from  the  surface,  S  is  the  molecular  speed  ratio 
u/fz  R  T*>  ,  is  the  temperature  of  the  incident  stream  and 

<3  is  the  angle  between  (J  and  an  element  of  surface  (the 
local  angle  of  attack).  The  wall  temperature  ratio  T* / Ho 
must  be  determined  from  energy  transfer  considerations  at 
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the  surface  and  will  be  discussed  later.  The  molecular 
speed  ratio  S  acts  like  a  Mach  number  and  is  in  fact 
related  to  the  Mach  number  through  5  -  • 

The  resulting  similarity  parameters  for  pressure 
distribution  about  a  body  in  free  molecule  flow  are  seen 
to  be 


S  s/a/  fi  y  Ik.  j  P  , 

!<*>  L 


(2.22) 


The  parameter  -f  depends  primarily  on  the  surface  material 
and  condition  and  generally  varies  between  .8  and  1.0.  The 
unsteady  condition  enters  through  the  motion  of  the 

boundary,  and  defines  the  reference  time  i0  .  The  Reynolds 
number  plays  no  role  other  than  that  of  serving  to  define  the 
flow  regime  according  to  the  Tsien  criterion,  //?frf0  >  IO- 

In  high  speed  free  molecule  flow  around  the  nose 
of  a  blunt  body  (  5  »  1  ,  s/v/3  of  order  of  1,)  Eq.  (2.21) 
can  be  reduced  to 


=•  5/  ^  ft 

■fs 

where 


(2.23) 


Thus  for  high  speed  free  molecule  flow  the  pressure  distri¬ 
bution  about  a  blunt  body  tends  to  become  independent  of 
s  (or  Mco)  ,  ~&l %,  and  f  .  One  need  only  simulate  here  the 
same  geometry  ( ft  )  as  in  the  continuum  hypersonic  flow  case. 


*Note  that  slightly  differing  values  of  P  can  readily  be 
accounted  for  in 
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This  independence  of  S  at  large  values  of  S  is  also 
shown  to  exist  theoretically  and  experimentally  for  the 
lift  and  drag  coefficients  about  various  bodies  in  Refs.  13 
and  14 . 

For  affinely  related  slender  bodies  of  differing 
thickness  ratios  T  >  the  form  of  Eqs .  (2.10)  and  ^2.22) 
leads  to  the  similarity  parameters 


Sr  , 


~TZj 

Too 


U± o 

L 


(2.24) 


At  very  high  molecular  speed  ratios  S  such 
that  5t»  1  even  though  T4l  ,  Eq.  (2.21)  reduces  to 


-f 

L  ulrz(z-o 


+  3 


ax  \ 
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u, 
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(2.25) 


Leaving  only  the  parameters 


s 


^  o 

Lr 


(/"to 

TT 


(2.26) 


as  the  required  similarity  for  S  T  >7  •  In  Ec[’  (2.25), 

w  represents  the  deflection  of  the  midplane. 

The  above  behavior  at  high  5  for  blunt  bodies 
and  at  high  ST  for  slender  bodies  is  exactly  analogous  to 
that  of  continuum  hypersonic  flow  as  discussed  in  Section  2.2. 
Thus  the  Mach  Number  Independence  Principle  appears  operative 
in  the  free  molecule  flow  regime  as  well  as  in  the  continuum 
hypersonic  region. 
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The  aerodynamic  flow  over  a  cool  body  in  high 
speed  flow  causes  an  aerodynamic  heat  input  rate  to 

flow  into  the  body.  This  ,  resulting  from  the  viscous 

and  boundary  layer  effects,  causes  a  subsequent  rise  in  the 
body  temperature.  As  is  shown  later  by  an  examination  of 
the  gas  body  boundary  condition  in  Section  3.1,  similarity 
of  the  subsequent  body  temperatures  requires  simulation  of 
the  Biot  number  parameter. 


T*  L  (2.27) 

KoTo 


w^erG  ft °  represents  the  reference  heat  conduction  coeffi¬ 
cient  of  the  body.  To  represents  a  reference  temperature 
(can  be  taken  equal  to  Ts  if  desired),  and  fA  is  the 
aerodynamic  heat  transfer  rate  over  the  body  surface. 

Since  -  A  (dT/d»)SA$  at  the  surface,  the  above  condition 

of  Eq.  (2.27)  can  be  reduced  to  a  condition  of 


(2.28) 


provided  the  non-dimensional  temperature  distribution  T 
is  the  same  at  corresponding  points  in  the  flow  field  close 
to  the  body  (boundary  layer).  This  will  occur  if  the  complete 
similarity  conditions  of  Eq.  (2.2)  for  supersonic  or  Eq.  (2.5) 
for  hypersonic  flow  are  fulfilled. 

To  search  for  less  restrictive  similarity  condi¬ 
tions  pertaining  specifically  to  the  heating  of  the  body, 
various  expressions  for  the  aerodynamic  heat  transfer  rate 
are  examined.  These  can  then  be  placed  into  Eq.  (2.27)  to 
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obtain  specific  combined  parameters  which  may  be  less 
restrictive  than  the  parameters  of  Eq.  (2.2)  or  (2.5)  plus 
Eq.  (2.28). 

a)  Stagnation  Point 

Vhe  most  severe  heating  rates  of  a  hypersonic 
vehicle  will  probably  be  experienced  at  the  stagnation  point 
of  the  body  or  wing  leading  edge.  The  high  stagnation 
temperatures  behind  the  shock  wave  of  a  blunt  nosed  surface 
in  hypersonic  flow  will  dissociate  the  air  which  may  sub¬ 
sequently  recombine  in  the  gas  phase  of  the  flow  or  at  the 
body  surface  if  the  surface  is  catalytic.  This  recombina¬ 
tion  process  may  appreciably  add  to  the  heat  transferred  by 
normal  molecular  conduction  under  certain  circumstances. 

The  boundary  layer  at  the  stagnation  point  will 
most  likely  be  laminar.  The  stagnation  point  heat  transfer 
rate  for  such  a  dissociating  gas  in  hypersonic  flow  is 
given  by  Fay  and  Riddell  in  Ref.  15  as 


(2.29) 


where  R  is  the  nose  radius  of  curvature,  is  the  amount 
of  the  stagnation  enthalpy  J, s  outside  the  boundary  layer 
involved  in  dissociation  and  Lc  is  the  Lewis  number  fD*/*. 

The  Lewis  number  for  air  does  not  vary  appreciably  and  is  in 
the  range  from  1.0  to  1.5. 

The  effects  of  dissociation  are  accounted  for 
principally  by  the  bracketed  term  of  Eq.  (2.29).  If  the 
dissociation  species  recombine  rapidly  in  the  gas  phase 
(equilibrium  flow)  n  =  .52,  and  if  recombination  proceeds 
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slowly  (frozen  flow)  n  =  .63.  If  the  flow  is  frozen,  the 
dissociated  air  may  diffuse  through  the  boundary  layer  and 
recombine  at  the  body  surface  if  the  surface  is  catalytic. 
This  latter  frozen-catalytic  wall  case  produces  almost  the 
same  heat  transfer  as  does  the  equilibrium  flow  case.  Only 
in  the  frozen-non-catalytic  wall  case  is  there  any  appre¬ 
ciable  reduction  in  .  In  such  an  event,  consideration 

of  such  effects  as  the  recombination  rates  in  the  gas  and 
at  the  wall  is  necessary  for  detailed  simulation  of  the 
heat  transfer  process.  These  effects  are  described  further 
in  Ref.  15  and  16. 

Making  use  of  the  hypersonic  approximations  for 
the  stagnation  point  quantities,  Eq.  (2.4),  and  placing  <%A 
of  Eq.  (2.29)  into  the  heat  conduction  similarity  parameter 
fAL./KoTo  will  result  in 


7aL 

K.7T 


?*UL  Pys 

4 3  To 


(2.30) 


where  L  is  now  a  characteristic  length  pertaining  to  the  nose 
region,  and  Jis  is  the  stagnation  enthalpy  of  the  flow 

if  ic  results  purely  from  velocity  conversion  into 
enthalpy) .  From  the  above  the  similarity  parameters  for 
aerodynamic  heating  become 


(2.31) 
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The  first  parameter  is  seen  to  be  a  combined  grouping  of 
the  individual  parameters  for  similarity,  that  is  Eq.  (2.5; 
plus  Eq.  (2.28),  and  as  such,  considerable  adjustment  can 
be  made  of  them.  The  brace  (  5  quantity  is  a  weakly 

varying  function  of  temperature  and  for  air  does  not  vary 
more  than  a  few  percent  over  the  range  to  be  encountered 
in  re-entry  heating  problems.  The  second  parameter  requires 
that  similar  wall  temperatures  be  maintained.  It  is  to  be 
noticed  that  for  initially  cool  bodies  in  high  stagnation 
temperature  air  streams  where  >  bhe  term  may 

be  neglected  in  comparison  with  Jt  ?  in  Eq.  (2.29).  This 
then  causes  the  second  parameter,  / Ts  ,  to  drop  out, 

leaving  additional  freedom  in  the  choice  of  To  for  the 
first  parameter  (provided  both  model  and  prototype  operate 
in  the  X,  cc  Ts  region)  . 

For  engineering  purposes  the  following  simple 
relation  has  been  suggested  in  Ref.  17  for  stagnation  point 
heating  in  a  hypersonic,  highly  cooled  boundary  layer 
(  Tu  <<  Xs  ),  in  equilibrium,  laminar  flow, 

/,  (2.32) 

where  Ct  is  some  constant  depending  on  the  gas.  This  is 
obtainable  from  Eq.  (2.29)  by  making  a  number  of  simplifying 
assumptions  and  seems  to  give  reasonable  estimates.  (This 
expression  also  applies  for  frozen  flow  provided  the  surface 
is  highly  catalytic.)  Placing  this  into  the  heat  conduction 
similarity  parameter  iaL/k.oT0  results  in  the  single  simil¬ 
arity  parameter  for  aerodynamic  heating. 
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c.JZZ  u* 


or 


(2.33) 


&  To 


The  second  form  of  this  is  seen  to  be  very  similar  to  the 
more  precise  term  of  Eq.  (2.31)  with  its  brace  |  ]  and 

bracket  [  J  terms  removed,  if  the  assumptions  that  f^s  "'/tT 
and  /z  are  made.  It  is  considerably  simpler  to  use 

Eq.  (2.33)  for  rough  effects  of  similarity. 

It  has  been  shown  in  Ref.  18  that  the  heat  trans¬ 
fer  distribution  /('f/>\rA<;  i-n  t^ie  vicinity  of  the  stagnation 
point  for  laminar  flow  is  relatively  insensitive  to  varia¬ 
tions  in  Mach  number  and  Reynolds  number,  but  rather  depends 
mainly  on  the  geometry.  Thus  for  engineering  purposes, 
simulation  of  the  laminar  heat  transfer  rate  around  a  blunt 
body  up  t:o  angles  of  about  80°  from  the  stagnation  point  is 
achieved  if  the  same  geometry  is  maintained.  This  however 
depends  on  the  flow  being  maintained  laminar  over  this 
extent.  The  question  of  whether  transition  to  turbulent 
flow  will  occur  in  this  extent  of  80°  is  intimately  tied  in 
with  the  Reynolds  number  of  the  flow.  There  seems  to  be 
some  evidence  that  for  many  flight  missions  of  the  prototype, 
laminar  flow  will  prevail  because  of  the  low  densities  of 
the  higher  atmosphere  and  the  highly  cooled  walls.  If  this 
cannot  be  assured,  it  seems  necessary  to  either  simulate 
the  Reynolds  number  individually  for  both  model  and  proto¬ 
type  as  well  as  some  other  quantities  such  as  roughness,  or 
else  attempt  to  artificially  induce  transition  by  trip  wires 
on  the  model. 

b)  Regions  Away  from  Stagnation  Point 

On  flat  plate-like  surfaces  and  at  regions  away 
from  the  stagnation  point  where  pressure  gradient  effects 
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are  insignificant,  the  following  flat  plate  expressions 
for  the  heating  rate  f A  are  often  used 


-  .331  Oti  & 
6  ~x 


Ur  -X> 


(laminar) 


(2.34a) 


h  t 

fA  =  ,03o(py)  (RC)  A-  (Ay  -Jw)  (turbulent)  (2.34b) 

* 


where  Ay  is  the  recovery  enthalpy  (  Ar  =J,y  +a,u/z)  •  These 
expressions  have  successfully  correlated  data  for  un — 
dissociated  flow  at  high  Mach  numbers  utilizing  local 
values  of  the  gas  properties  just  outside  the  boundary 
layer.  They  are  also  applicable  to  the  region  around  a 
blunt  body  aft  of  the  stagnation  point  provided  the  pressure 
gradient  does  not  vary  too  rapidly, (Ref .  19). 


Making  use  of  the  hypersonic  approximations  and 


placing  the  above  heating  rates 

into  the  heat  conduction 

similarity  parameter 

L  /K0  ~fo 

results  in  the  parameters 

'h 

-A.  a;  (p*~)  -/v 

IS  -T- 

Tw 

> 

(laminar)  (2.35a) 

Ko  yCr  To  Ts 


•  *  7* 

A  (fc)  (Pv) 

(turbulent) 

(2.35b) 

/Co  "7*? 

7> 

The  above  combined  similarity  parameters  are  quite  similar 
in  form  to  those  derived  for  the  laminar  stagnation  point 
heating  rate,  Eq.  (2.31).  Again,  for  highly  cooled  walls. 
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the  second  parameter  ~^>/~rs  can  be  neglected  in  the  similar¬ 
ity.  It  is  to  be  noted  that  a  small  conflict  will  result 
if  the  above  combined  parameters  are  attempted  to  simulate 
simultaneously  both  the  laminar  and  turbulent  regions  of 
flow  because  of  the  different  powers  of  Pc  •  However 
the  discrepancy  is  not  great  and  some  inaccuracy  could 
probably  be  tolerated  without  having  to  resort  to  matching 
,  &/ 'Ko  >  and  P-^  separately.  (Of  more  importance 

would  be  insuring  that  transition  occured  at  the  same  point 
on  model  and  prototype.) 

c)  Sharp-Nosed  Slender  Bodies 

For  a  family  of  affinely  related  slender  bodies 
of  thickness  ratio  'P  ,  Hayes  and  Probstein  in  Ref.  2  have 
shown  through  an  analysis  of  the  hypersonic  laminar  boundary 
layer  that  the  laminar  heating  rate  can  be  expressed  as. 


h  =  Ml  A  (%U;  a U,  *  ,  r,  rr,  £- 

ifTur  Lr 

V  Ps 

where  A  is  a  function  of  the  parameters  shown,  and  the 
stagnation  enthalpy  J\  r  has  been  taken  as  J,f  ~  Ux/z  . 
Placing  the  above  into  the  heat  conduction  similarity  para¬ 
meter  fa  L/KoTo  results  in  the  combined  parameters. 


Uh 

L 


■) 


(2.36) 


A. /p+ul  ruL  ,  7^  ;  MaoT  t  y 

Ko  V  'C-fs  To  Is 


Av- 
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L 


(2.37) 
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It  is  seen  that  the  first  parameter  is  a  more  relaxed  version 
of  the  previous  flat  plate  parameter,  Eq.  (2.35a),  since 
here  'Y  also  appears  in  the  combination.  (The  difference  in 
the  /V  power  is  not  significant  since  Pr  is  required 
to  be  maintained  separately  anyway  in  Eq.  2.37).  The 
second  wall  temperature  parameters  are  identical, 

and  the  remaining  parameters  of  Eq.  (2.37)  are  already 
required  to  insure  that  aerodynamic  pressure  similarity  is 
being  maintained  over  these  slender  bodies. 

As  mentioned  in  Section  2.2b,  viscous  interaction 
between  the  boundary  layer  and  the  flow  field  may  play  a 
significant  role  at  these  hypersonic  speeds.  To  maintain 
the  proper  viscous  interaction,  one  would  need  to  add,  as 
in  the  pressure  distribution  case,  the  interaction  parameter. 


T  =  (2.38) 

to  the  parameters  given  by  Eq.  (2.37).  For  very  high 
Mach  numbers  where  A?*®  "^  »  i.  ,  the  viscous  interaction 
parameter  of  Eq.  (2.38)  is  replaced  by 


/ 

r^Ul 


(2.39) 


and  the  Wo*  Y  parameter  drops  out  of  the  grouping  given 
by  Eq.  (2.37). 

For  heat  transfer  effects  generally,  the  real 
gas  effects  are  expected  to  play  a  larger  role  than  for  the 
pressure  distribution  effects.  It  would  then  probably  be 
wise  to  see  that  at  least  identical  temperatures  are 
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maintained  on  the  model  and  prototype.  For  this  case,  the 
parameters  become  those  of  Eq.  (2.16)  plus  the  combined 
parameter. 


As  I  f>+UL  Pr$  r  HL 

Afo  *  /*s 


(2.40) 


If  further,  the  same  pressures  as  well  as  the  same  temper¬ 
atures  are  to  be  maintained  in  order  to  insure  the  same 
equation  of  state,  then  similarity  would  require  the  para¬ 
meters  of  Eq.  (2.17)  plus  Eq.  (2.40).  Even  then,  it  still 
might  be  necessary  to  consider  recombination  effects  for 
complete  similarity. 

To  account  for  the  presence  of  blunt  nosed  slender 
bodies,  the  parameters  of  Eq.  (2.18),  Mj  A  or  At*  A  ^  } 

should  be  added  to  those  already  existing. 

d)  Free  Molecule  Flow  Regime 

Utilizing  concepts  from  the  kinetic  theory  of 
gases,  the  heat  transfer  rate  in  free  molecule  flow  can  be 
determined  analogously  to  the  pressure  distribution.  Follow¬ 
ing  Ref.  20  this  heat  transfer  can  be  expressed  as. 


f/i  = 

vjw  S3 


X-  -  XX  Xl 

L  Ur-0  x. 


f  JlV  S  5/^/3  [l  +■  (  S  S/v  •)])  -  ± 


-  f  S/n/l/3 


-  5  s/W 


(2.41) 


where  as  before,  S  is  the  molecular  speed  ratio,  y  is 
the  specific  heat  ratio  and  is  an  accomodation  coeffi¬ 

cient  depending  on  the  surface  and  is  usually  in  the  neighbor¬ 
hood  of  unity. 


28 


Placing  the  above  into  the  heat  conduction 
similarity  parameter  'fAL/Ko'fe  results  in  the  following 
general  similarity  parameters. 


«e  p^'A  ,  stmt  ,ik  ,r,s,  uh.  <2.«) 

■sV.rr  1  L 

The  first  parameter  is  the  basic  heat  conduction  requirement. 

It  is  interesting  to  note  that  it  can  be  rearranged  to  give, 

ZZe  (P<t°UL\  Prco  Aco  JL-  (2.43) 

V  A**  /  5  3  Ko  X^% 

thus  resembling  strongly  the  hypersonic  continuum  parameters 
of  Eq&  (2.31),  (2.35),  and  (2.37).  Note  that  since  S  is 
required  similar  separately,  it  can  be  dropped  from  this 
first  parameter.  The  second  and  third  parameters  of  Eq.  (2.42) 
are  already  required  for  simulation  of  the  pressure  distri¬ 
bution,  (see  Eq.  2.22),  while  the  J  requirement  is  a 
relatively  weak  one  mainly  describing  whether  a  monatomic 
or  diatomic  gas  is  present.  The  final  requirement  is  that 
the  speed  ratio  S  be  maintained  separately  as  well  as  in 
combination  with  S'  s/r,  0  . 

In  high  speed  free  molecule  flow  around  the  nose 
of  a  cool  blunt  body  (  S  >>  ^  >  s/n/3 of  order  JL  }  5  /3  )> 

Eq.  (2.41)  can  be  reduced  to 


which  when  placed  into  the  heat  conduction  parameter  ~0 

yields  the  single  heating  similarity  requirement , 


K.  To 


(2.45) 


The  sin/3  term  has  been  dropped  in  the  above,  since  we  are 
considering  bodies  of  similar  geometry  (  (i  distribution) 
here.  Thus,  it  is  seen  that  the  heating  characteristics 
tend  to  become  independent  of  S  ,  7*/voP  and  y  for  this  high 
speed  cool  blunt  body  case. 

For  affinely  related  cool  slender  bodies  of 
differing  thickness  ratios  T  in  moderately  high  speed  flow 
(  5  >3  j  «  s /3  ),  Eq.  (2.41)  can  be  reduced  to 


‘tJW  S 


f  _  (5 smfi) 

j  e  -h\ffT  Sf/t/fi 


(2-46) 


The  form  of  Eq.  (2.46)  then  leads  to  the  similarity  parameters 


fa,  U  L  T  ,  $r  >  $  (2.47) 

KoTo  LT  '  ~ 

where  the  first  parameter  now  has  the  thickness  ratio  T 
combined  with  it.  The  5r  is  the  familiar  hypersonic 
similarity  parameter  encountered  in  dealing  with  pressure 
distributions,  while  the  next  two  parameters  insure  simil¬ 
arity  of  thickness  distribution  (family),  and  subsequent 
deflect  ions . 

At  very  high  molecular  speed  ratios  5  ,  such 
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that  S 'Y  »  3-  even  though  'Y  l.  1  ,  Eq.  (2.46)  reduces 

to  Eq.  (2.44)  which  then  leads  to  the  dropping  of  the 
parameter  9T  from  eq.  (2.47)  leaving  only 

^ef)CoUlL'T  ,  $  ,  ,  Uio_  (2.48) 

~~fUTo  LT  L 

Thus  the  heating  behavior  parallels  the  pressure 
distribution  behavior  discussed  in  Section  2. 2d,  and  also 
exhibits  the  characteristics  of  the  Mach  Number  Independence 
Principle . 

In  dealing  with  heating  in  free  molecule  flow, 
the  above  discussed  aerodynamic  heating  rates  are  only  one 
of  the  inputs  to  the  system.  Since  these  are  in  general 
quite  small,  other  heating  inputs  of  comparable  or  even 
greater  magnitudes  must  also  be  considered  when  dealing  with 
heat  transfer  in  free  molecule  flow.  Three  of  these  that 
assume  major  significance  are  the  solar  heating  rate  %sol 
the  earth  heating  rate  ^eav  >  and  the  radiation  output 
from  the  body  wall  itself,  C*,  r  .  These  inputs  and 

outputs  bring  in  the  additional  heating  rate  parameters 

fs°L  L  f  f ear  L  >  €o  To  *L  (2.49) 

K°To  K»To  Ko 

into  the  problem.  In  their  simulation,  the  proper  distri¬ 
bution  of  these  over  the  surface  of  the  body  is,  of  course, 
required.  Often  times  the  above  parameters  may  be  of  much 
greater  importance  than  the  free  molecule  flow  aerodynamic 
heating  parameters  previously  discussed. 
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3. 


HEAT  CONDUCTION  SIMILARITY 


3.1  Non-ablating  Body 

The  equations  governing  the  distribution  of  heat 
within  a  body  are 


These  equations  are  similar  to  the  equations  treated  in 
Ref.  1  except  that  now  the  heat  flux  from  the  boundary 
layer  is  written  as  ^  in  place  of  (A  d7!?}) I  .  Eq.  (3.1) 
governs  the  conduction  of  heat  within  the  body  subject  to 
the  boundary  condition  of  Eq.  (3.2)  which  expresses  the 
balance  of  heat  flux  at  the  heated  surface.  As  in  Ref.  1 
an  insulated  wall  boundary  condition  (  K^T/dt)  o  )  will 
not  introduce  any  additional  parameters  to  the  equations 
above . 

Non-dimensionalizing  by  introducing  the  following 
quantities : 


X  = 

X 

L 

>  t 

-  r  j 

i  -  X 

L 

T  = 

T 

To 

,  * 

-Ms* 

II 

(3.3) 

= 

K_ 

Ko 

j  C?s 

- 

cTo  J 

€  */  = 
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the  above  equations  become 


where  TC„  —  Ko/?g  .  The  resulting  non-dimensional 

parameters  required  for  similarity  in  the  temperature  dis¬ 
tribution  T  are 


frL  }  €oH  L  )  JToio  j  Kg 

~KoTc  K O  O- 


> 


(3.6) 


The  above  parameters  are  similar  to  the  parameters 
derived  in  Section  2.2  of  Ref.  1,  where  the  parameter 
f*L  /k>t0  ,  which  may  be  identified  as  the  Biot  number, 

corresponds  to  the  A^/Ko  parameter  there.  It  should  be 
noted  however  that  replacing  faL/KeTo  by  A^/Ko  requires 
the  simulation  of  the  temperature  gradients  in  the  boundary 
layer  and  these  gradients  generally  depend  on  reproducing 
Moo  and  separately.  The  parameter  also 

depends  on  /^oo  and  through  ,  but  as  indicated  in 

the  preceding  section  W*,  and  f?eao  appear  in  combined 
forms  and  in  many  cases  (for  sufficiently  high  )  the 

dependence  on  may  be  dropped.  The  actual  form  that 

the  Biot  number  parameter  ^-fiL/KoH  does  take  on  for  various 
aerodynamic  theories  has  been  discussed  in  Section  2.3.  It 
is  to  be  noted  that  the  heat  flux  ^  in  the  Biot  number 
parameter  may  be  of  either  aerodynamic  or  non- aerodynamic 
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origin. 


because  of  the  resulting  higher  surface  temperatures  antici¬ 
pated  in  hypersonic  flight,  is  generally  of  considerable 
importance.  If  the  same  materials  and  temperatures  are  to 
be  utilized  reproduction  of  this  parameter  is  possible  only 
if  the  length  scale  ratio  LP  /l*  =  4-  .  Utilization  of 

the  same  materials  and  temperatures  is  mandatory  not  only 
for  simulation  of  the  heat  transfer  problem  (matching  Kg 

c  )  but  also  for  structural  deformation  behavior  as 
will  be  seen  in  the  next  section.  The  radiation  parameter 
thus  imposes  a  stringent  restriction  on  the  length  ratio 
Lr/L»i  and  presents  a  basic  obstacle. 

Several  possibilities  exist  in  dealing  with  this 
obstacle.  In  one  it  may  be  possible  to  control  the  emmis- 
sivity  by  such  means  as  coating  the  surface  of  the  structure 
and  thus  not  restricting  the  value  of  £  by  a  choice  of 
materials.  In  this  case,  for  the  same  materials  and  temper¬ 
atures,  the  scale  ratio  as  determined  from  £0Ti  L/fo  will  be 


Lp_ 

L-  A7 


(Vo)p 


(3.7) 


As  (fo)/YI  can  never  be  greater  than  1,  there  is  a  limit 
to  the  scale  ratio  Lp /l*  that  can  be  satisfied  by  this 
means,  particularly  if  the  prototype  itself  is  operating  at  high 
emissivities  (£^)F  .  In  utilizing  a  surface  coating  to 

obtain  the  proper  £c  ,  it  should  also  be  determined  whether 
this  coating  affects  the  heat  flux  . 

A  second  means  of  alleviating  the  radiation 
obstacle  occurs  when  considering  flight  near  the  equilibrium 
temperature.  The  equilibrium  temperature  is  achieved  after 
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a  sufficiently  long  period  of  heating  when  the  heat  flux 
input  balances  the  radiative  heat  loss.  Transfer  of  heat 
by  conduction  is  negligible  (i.e.,  (K  ^rhi)g-  o  in  Eq.  3.2), 
and  the  parameters  of  Eq.  (3.6)  are  replaced  by  the  following 
parameters 


U  >  (3.8) 

<f07TV 

Therefore  in  this  case  no  condition  on  the  length  scale 
appears  directly.  The  parameter  %  n  / is  particularly 
applicable  to  a  hypersonic  glide  vehicle  which  by  nature 
of  its  long  flight  time  relies  principally  upon  radiation 
as  a  cooling  mechanism. 

A  further  means  of  dealing  with  the  radiation 
obstacle  may  be  seen  by  rewriting  Eq.  (3.5)  in  the  following 
non-dimensional  form 

(<ell)  =  =  fcflL  (3.9) 

and  thus  directing  attention  to  the  net  heat  flux  £/w=r 
into  the  structure  rather  than  treating  fA  and  the  radiation 
loss  Cw<r  tJ*  as  two  separate  quantities.  The  parameters  of 
Eq.  (3.6)  would  then  be  replaced  by 

f*£rL  7V. -to  Kg  j  Crg  (3.10) 

TTtT  f 

This  procedure  assumes,  however,  that  and  Xj  are 

known  in  advance  (see  further  the  discussion  in  Section  5.1e). 
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Having  estimated  the  required  'J.net  for  a  particular 
test,  it  may  be  applied  by  such  means  as  radiant  heaters 
or  a  heated  air  stream.  Use  of  radiant  heaters  may  provide 
similarity  in  heat  transfer  and  allow  large  scale  ratios 
but  will  not  account  for  any  air  loads.  A  heated  air  stream 
at  a  low  may  provide  similarity  in  heat  transfer  but 

may  incorrectly  scale  the  air  load  distribution.  Tests 
of  this  second  type  have  been  conducted  on  multiweb  wings 
by  the  NASA  (see  Ref.  21,  for  example). 

Heat  transfer  similarity  has  been  discussed  above 
for  a  general  body  (Eqs.  3.1  and  3.2).  Some  less  restrictive 
relationships  can  be  obtained  however  when  considering  specific 
configurations  such  as  thin  skins  and  built  up  structures. 

The  general  ideas  involved  here  are  discussed  in  Ref.  1. 

For  example,  considering  a  thin  plate  of  thickness  which 
has  a  very  small  temperature  variation  through  the  plate  in 
the  h  direction,  the  governing  heat  conduction  equation 
becomes. 


=  f.c' 


(3.11) 


Non-dimensionalizing  this,  the  appropriate  similarity 
parameters  are, 

tl£ 1  ,  ±2LA  ,  xjdi.  s  0.12) 

Ko  To  So  K  o  2  o  L- 

The  above  are  somewhat  less  restrictive  than  the  parameters 
of  Eq.  (3.6)  since  they  allow  the  thickness  ratio  'X ~  So/L 
of  the  plate  to  be  varied  in  addition  to  the  original 
quantities.  However,  both  the  prototype  and  model  plates 
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must  behave  as  thin  plates  with  little  temperature  variation 
in  the  Z  direction. 

Even  further  reduction  in  the  similarity  parameters 
can  be  achieved  if  it  is  assumed  as  well  that  the  heat 
conduction  in  the  x  and  'f'  direction  is  also  very  small. 

In  that  case  the  parameters  of  Eq.  (3.12)  reduce  still 
further  to 


f'A'ko  j  Go  ^ 0  ,  C -p  ,  £ #  )  $ 

Pg  CK  SoTo  Pg  CTo  So 


Such  an  assumption  has  been  made  for  the  heat  conduction 
of  thin  two  dimensional  wings  in  Ref.  22. 


3.2  Ablating  Bod^ 


When  the  body  temperature  is  sufficiently  high 
the  material  may  undergo  a  change  of  phase  by  either  melting, 
accompanied  with  possible  vaporization,  or  by  sublimating. 
The  process  of  absorbing  heat  by  removal  of  the  surface 
material  is  known  as  ablation  and  is  of  considerable  current 
interest.  It  is  likely  to  occur  at  the  noses  or  leading 
edges  of  high  speed  vehicles  where  the  aerodynamic  heat 
transfer  rates  are  apt  to  be  very  large. 


The  ablation  process  of  a  melting-vaporizing 
material  near  a  stagnation  point  may  be  treated  by  studying 
the  laminar  boundary  layer  equations  for  a  melt  layer, 
which  have  the  following  forms  (see  Refs.  23,  24,  and  25) 


dfUsi-  )  P  ( <va  )  —  o 

d  x 


(3.14) 
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1 Z.  * 

9 X  ^2  1 


=  O  (3. 16) 
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fud'T  4-sv- 3 T  1 
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l  d7  32  J 
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d *  1 

where  j  =  0  for  a  two-dimensional  flow  and  j  =  1  for  a  body 
of  revolution. 

The  inertia  forces  are  small  compared  to  the 
viscous  forces  and  hence  are  usually  neglected  above  (/ ??t 
Also  the  viscous  dissipation  terms  in  the  energy  equation 
can  be  neglected.  The  accompanying  governing  equation  for 
heat  conduction  in  the  body  (Eq.  3.1)  is 


where  the  surface  is  a  distance  s  from  the  initial  pre¬ 
ablation  position  (see  Fig.  3).  The  boundary  conditions 
are 

at  z  =  h  ( t )  : 

f/f  C^s-^)  ~  ~  ^  fj \  +  -h  £-<r  / 

(3.19) 

T4  =  ■  ¥1  >  i  r=  % 

.j —  -dA  r  u  lA  — 

Pi  **  Pl. 
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at  z  =  s(t)  : 


where  sr*L  and  ^ /v  are  the  mass  rates  of  loss  due  to  melting 
and  vaporization,  respectively,  L  r~\  and  L^r  are  the  heats 
of  melting  and  vaporization  respectively,  is  an  empirical 

factor  associated  with  the  heat  blockage  due  to  vapor  injec¬ 
tion,  7^  and  Zy  are  the  melting  and  vaporization  temper¬ 
atures,  respectively,  is  the  stagnation  enthalpy,  and 

4  is  the  ablation  velocity.  The  aerodynamic  heat 

transfer  fA  ,  shear  T*  and  pressure  -f*  are  assumed  known 
quantities  for  a  non-ablating  material,  evaluated  at  the 
gas-melt  interface,  and  were  discussed  in  Section  2.  (The 
shear  stress  is  directly  related  to  the  heat  transfer 

through  Reynold's  analogy,  X  /t*  =  (rYsf*«  /(Js-Ai) 
where  Pys  is  the  stagnation  Prandtl  number  of  the  gas. 

Non-dimensionalizing  by  introducing  the  following. 


u  = 


u_ 

Vi 


jr  - 

Vo 


, 


t- 

% 


i 


~T  - 


X 

To 


■L  — 


t 
£  0 


(3.21) 


^  Po  ‘  Cfi 


p  ■=  Pr 

r*  %  J 


=  CfL 
Cfo 


Cw 
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equations  (3. 14)  to  (3. 18)  become 


(3.22) 


(3.23) 


££  ^  o 


(3.24) 


) fe  fr  ’  k  &  U ')  *h (*’ r)  <3’26) 

where  7T.  =  Ko/?0cr.  is  the  thermal  diffusivity.  Since  stag¬ 
nation  point  ablation  is  being  considered,  the  characteristic 
length  L  in  this  case  is  appropriately  the  nose  radius  R. 
The  reference  quantities  Ko ,  f>0  )  C-p0  refer  to  reference 

properties  on  the  body. 
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are 


The  boundary  conditions  in  non-dimensional  form 


at  z  -  R(t ) : 

^  fA  L  J  _  Cf0  L-at  +  ^  ^ 7c  L)  £  f  ~Ti 

“  '  (3.27 ) 

u  3 u  'Ya  L  -p  -  ^  -  "7£ 

'  7  ?  /Wc  T'a  J  7^ 

-  /_kL 1 +  u  zi  -fjL  ){?l) /-rk ) 

^  =  '  X'.-io/WoL)  d  i  ??  ITT Ao/(Vdl/\  1A 

at  z  5  ( t )  : 


„  _  (3.28) 

-KcfiZ)  =  '  XfifZZ)  +■  Anu-QlL)  ( L”  \ 


The  important  unknown  non-dimensional  quantities 
to  be  solved  for  from  Eqs.  (3.22)  to  (3.28)  above,  are  the 
ablation  rate,  the  fraction  of  melt  vaporized,  and  the 
temperature  distribution  within  the  body 

mlCtU.  ,  ^  ,  Tg  (3.29) 

Ko 

The  parameters  on  which  these  unknowns  depend  are 
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obtained  from  Eqs.  (3.22)  to  (3.28)  as  follows, 

fAL  }  eoTZL  ;  TT.io  ,  ,  M-  ,  T4  L 

Ko7o  K*  L2~  "To  C-fo  71  j^eVo 

foL_  }  VoL  J  7v  y  <7  (3.30) 

/'•Vi  7T.  *  To  K 

KB  ,  ?-p0  ,  Pg  ,  *L-  >  j  ?L.  >  fi-L 

The  above  are  applied  to  high  speed  stagnation  point  flow  by 
setting  the  reference  pressure  -f0  equal  to  the  stagnation 
pressure,  then  substituting  f or from  Eq.  (2.6)  and  for 
from  Reynolds  analogy 


-fo-=1>s  -  fcV1- 

fA  *  f*>  u'  *'”2/3  (3.31) 

V) 

TA  =  Cz  fa  (f$)_ 

V  R 

where  Cz  is  some  numerical  constant. 

The  similarity  parameters  of  Eq.  (3.30)  finally 
become,  after  slight  rearrangement. 


faL 

Ko  75  ; 


C0-To3L  ,  7u>^o_  } 

Ko  '  L *- 


Tv\  L  /V; 

~r°  C-fi  To 


Vi 

(f*s)  KoToL 
CJ  TCo 


,  Mi\ 

ftoKo  JCo 


T 

To 


) 


L*  +  n  (-4,  -X) 

c*  % 
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(i  ,  *.  ,  4.  ,  h  .  *T  4  ,  ft-  ,  p"-  (3.32) 

The  above  are  the  main  ablation  similarity  parameters.  The 
first  three  parameters  are  the  same  as  derived  previously 
in  Section  3.1  for  the  non-ablating  body.  The  following 
five  parameters  refer  to  the  melting  ablation  process. 

Taking  them  in  order,  the  fourth  parameter  serves  to  define 
the  reference  temperature  To  ,  the  fifth  gives  a  condition 
on  the  latent  heat  of  melting,  the  sixth  on  the  shear  forces 
acting  on  the  melted  layer,  the  seventh  on  the  pressure 
forces,  and  the  eighth  serves  to  define  the  reference  veloc¬ 
ity  Vo  associated  with  the  flow  of  the  molten  layer.  The 

next  two  parameters  refer  to  the  vaporizing  process  and 

•k 

give  conditions  on  the  vaporizing  temperature  and  effective 
heat  of  vaporization  respectively,  and  the  eleventh  is  a 
condition  on  geometry.  The  remaining  parameters  require 
similarity  in  variation  with  temperature  of  the  material 
properties  both  in  the  solid  and  the  liquid  range.  As 
given  in  Eq.  (3.32)  the  parameters  apply  to  a  combined 
melting  and  subsequent  vaporization  process.  If  only 
melting  is  present  (as  in  metals  or  at  early  times  before 
a  sufficient  vaporization  temperature  is  achieved),  the 
vaporization  parameters  drop  out.  If  the  ablation  is  a 
sublimation  process,  then  the  melting  parameters  drop  out. 
Also,  if  all  the  melt  is  removed  immediately  as  it  is  formed, 
the  liquid  layer  does  not  enter  the  problem,  and  the  6th,  7th, 
and  8th  parameters  above  which  relate  to  the  flow  of  this 


^Note  though  that  the  vaporization  temperature  of  some 
materials,  such  as  glasses,  may  themselves  be  complex  functions 
of  other  parameters  of  the  problem,  for  example,  pressure. 
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layer  would  not  be  present. 

Reproduction  of  the  parameters  of  Eq.  (3.32) 
will  simulate  a  detailed  study  of  the  ablation  process.  In 
particular,  the  non-dimensional  quantities, 

/mL  C-p,, L  }  /rrt</  t  7^  (3.33) 

Ko  /mu 

will  be  the  same  for  both  prototype  and  model  at  corres¬ 
ponding  non-dimensional  times  t.  In  many  instances,  the 
structural  designer  is  not  concerned  with  the  details  of 
the  ablation  process  but  rather  is  interested  in  simulating 
only  the  temperature  distribution  T  within  the  body,  and  in 
knowing  the  rate  of  ablation  ^  .  In  this  case,  the  net 

balance  of  heat  flux  at  the  surface  z  =  s  is  of  interest. 
With  the  use  of  Eqs.  (3.19)  and  (3.20)  this  net  heat  balance, 
replacing  Eq.  (3.2),  is 


}a  '  ^ 


~Af 


*  *7  (As-X) 


1 


(3.34) 


Assuming  steady- state  ablating  conditions  are  reached,  it 
can  be  shown  that  (see  Ref.  25,  for  example) 


(kl  a) 


—.60  /yvj^  ^  I*  ~  ”Ai ) 


(3.35) 


%  ^ 


(Tm-T*) 


(3.36) 
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Utilizing  Eq.  (3.35)  and  (3.36),  the  heat  flux  relation, 
Eq.  (3 .34)}  becomes  in  non-dimensional  form. 


^ lCp0L 


PlL\ 

<»To) 


\  lo  / 

[cij  (Tm'~Tq)+  L"*-  to  CfL  (%-rM)  +  n.(4s-A*)]f^ 

I  Oi  To  (  e«f%V 

K  9  A 


The  above  gives  the  non-dimensional  melting  rate  for 
steady  state  ablation.  It  represents  a  combined  grouping 
of  the  individual  ablation  parameters  given  in  Eq.  (3.32), 
and  hence  also  gives  some  indication  as  to  which  are  the 
more  significant  of  these  parameters.  If  a  judicious  guess 
for  the  unknown  sr*v  appearing  in  it  can  be  made,  the 

steady-state  ablation  rate  is  then  known.  Also,  the  rough 
temperature  distribution  within  the  body  could  be  obtained 
from  Eq.  (3.36)  as. 


Tm  on  the  surface 


slope 


_  JL  (~L  Cr.L) 
'  Ko  / 


at  the  surface 

(3.38) 


Tn  in  the  interior 


In  dimensional  form,  Eq.  (3.37)  is  often  written 


(3.39) 
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where  H  ^  is  defined  as  the  effective  heat  of  ablation 

G  jl  jC 

in  units  of  BTU/  lb,  and  is  given,  for  the  case  considered 
above,  by 


H, 


(Tn-Tc)  *  L„  +.toCK  (%-T„) 

j[  —  (T  ~Tv 

_____ 


(3.40) 


It  is  seen  that  radiation  may  increase  He££  considerably. 
Furthermore,  if  the  radiation  loss  is  neglibible  ( €w r %  »  O  ) , 

He££  becomes  independent  of  the  length  scale  L  which  appears 
through  the  term  ft\  .  Efficient  ablators  generally  operate 
with  considerable  vaporization  present,  and  so  the 

[  J  term  generally  predominates  over  the 

preceding  terms.  The  above  expression  for  Hg££  is  fre¬ 
quently  used  by  structural  designers  to  estimate  roughly 
ablation  shield  requirements  (see  Ref.  26) .  However,  it 
should  be  born  in  mind  that  Eqs .  (3.37)  and  (3.40)  represent 
steady-state  ablation  rates,  and  that  is  estimated. 

For  a  more  accurate  detailed  simulation  of  the  transient 
ablation  process,  the  parameters  of  Eqs.  (3.32)  should  be 
used . 
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4. 


STRESS  AND  DEFLECTION  SIMILARITY 


The  high  structural  temperatures  and  air  loads 
associated  with  flight  at  high  Mach  numbers  may  no  longer 
justify  the  assumption  of  an  elastic  structure  made  in 
Ref.  1,  and  additional  structural  and  material  considerations 
are  necessary.  Structural  deflections  may  extend  into  the 
plastic  range  in  which  case  material  yield  and  ultimate 
stresses  and  plastic  stress-strain  behavior  must  be  con¬ 
sidered.  Furthermore,  the  assumption  of  small  deflections 
may  have  to  be  removed  and  it  may  be  necessary  to  scale  the 
structural  deflections  in  the  same  scale  as  the  character¬ 
istic  dimensions  of  the  body.  At  high  temperatures  the 
creep  and  fatigue  characteristics  of  a  structure  may  have 
to  be  considered  which  introduce  time,  temperature  level 
and  stress  level  as  additional  parameters.  For  convenience, 
the  stress  and  deflection  problem  is  divided  into  two  major 
categories,  namely,  short-time  loading  and  long-time  loading. 

4.1  Short-Time  Loading 
a)  General 

This  category  is  concerned  with  loads  of  sufficiently 
short  duration  such  that  the  stress  deformation  behavior  of 
the  material  is  not  influenced  appreciably  by  the  time  under 
load,  i.e.,  creep  and  fatigue  effects  are  absent.  The 
material  behavior  can  then  be  described  solely  by  means  of 
stress-strain  curves  at  the  particular  operating  temperature. 

It  should  be  noted  that  for  very  short  duration  loads  (high 
strain  rates),  these  curves  tend  to  become  additionally 
dependent  on  strain  rate  £  .  The  stress  deflection  be¬ 

havior  considered  here  will  be  at  stress  levels  between  the 
yield  point  and  the  ultimate  stress.  The  parameters  for 
the  elastic  range  have  been  derived  previously  in  Ref.  1. 
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In  this  regime  nonrecover able  deformations  may 
result  which  may  alter  the  steady-state  deflections  of  the 
wing,  the  internal  structural  geometry,  and  the  wing  profile. 
These  deformations  may  subsequently  alter  the  air  load 
distribution  acting  on  the  structure.  Plastic  buckling 
may  also  occur  in  members  subjected  to  compressive  forces. 
Essentially  then,  this  category  is  an  extension  of  the 
elastic  structure  study  of  Ref.  1,  Section  2.3,  into  the 
plastic  range  and  the  additional  parameters  which  must  be 
included  relate  to  the  short-time  plastic  stress-strain 
behavior  of  a  material.  The  similarity  parameters  required 
for  the  same  stress  distribution  are  thus 

shape  of  stress-strain  curve, 

P'L1  .  PjzlL  ,  £  ,  “  ,  (4,1) 

/  2.  /  1 
00  to  &o 

The  first  three  parameters  above  relate  to 
static  stress-strain  behavior,  the  fourth  to  the  body 
temperature  distribution,  the  fifth  to  scaled  structural 
deflections,  the  sixth  to  structural  dynamic  effects,  and 
the  seventh  to  gravitational  forces.  The  last  three  are 
conditions  on  material  properties  which  are  automatically 
fulfilled  if  the  same  materials  and  temperatures  are  utilized. 
For  convenience,  the  reference  stress  r0  may  be  taken  as  the 
yield  stress  of  the  material.  Note  that  for  an  elastic 
body  the  condition  on  the  shape  of  the  stress -"Strain  curve 
above  is  automatically  met. 

The  simulation  of  the  general  shape  of  the  stress- 
strain  curve  is  a  very  difficult  task  for  anything  except 
the  same  material  at  the  same  temperature.  An  attempt  at 
this  can  be  made  though,  through  the  Ramberg- Osgood 


48 


approximation.  It  has  been  shown  that  for  many  metals  at 
various  temperatures  and  various  soaking  times  the  following 
Ramberg- Osgood  expression  is  a  reasonable  approximation  (Ref. 27) 

<4‘2) 

where  £7  is  the  stress  corresponding  to  a  secant  modulus 
of  . 7E  and  n  is  the  shape  factor  of  the  curve.  Both  E 
and  <T7  are  now  functions  of  temperature  T.  Non-dimen- 
sionalizing,  Eq.  (4.2)  becomes 


where 


(4.3) 


In  the  above,  (<r.7\  is  the  cC7  at  some  reference  temperature. 

Thus  the  entire  shape  of  the  stress-strain  curve  is 
reasonably  well  reproduced  if  the  following  parameters  are 
duplicated 


ft  j  £  ■  ft  ,  (T7  ,  1  ,  ^<-T  (4-4) 

E  o  (^7)e  ft 

The  first  and  second  conditions  above  were 
previously  derived  for  a  body  in  the  elastic  region;  the  re¬ 
maining  ones  are  now  additional  conditions  for  the  plastic 
portion  of  the  stress- strain  curve.  The  parameter  r°/(t£) 
serves  to  define  the  reference  stress  level  <T0  ,  while  the 

<£>  condition  requires  the  to  have  the  same  variation 

with  temperature  on  model  and  prototype.  The  parameter 
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is  introduced  to  account  for  the  extent  of  the  plastic 
region  and  may  be  dispensed  with  if  both  model  and  prototype 
operate  below  their  ultimate  strengths.  With  regard  to  the 
shape  parameter  n,  for  4340  steel,  stainless  steel  W  and 
Inconel  X,  for  example,  n  >  10  up  to  1200°F  and  n  ~  10 
approximates  many  of  the  aircraft  metal  alloys  at  room  and 
elevated  temperatures.  The  shape  parameter  n,  being 
approximately  the  same  for  many  aircraft  materials,  may 
therefore  be  dispensed  with.  Under  all  the  above  conditions, 
the  first  two  parameters  of  Eq.  (4.1)  may  then  be  replaced 
by  the  following: 


HL  >  1 L  )  7  (4.5) 

Eo  (t?)0 

where  it  has  been  assumed  that  all  stresses  are  acting 
below  the  ultimate  stress  of  the  material.  An  additional 
condition  is  necessary  relating  to  the  normal  stress  boundary 
condition  on  the  structure.  At  the  surface  of  the  body,  the 
normal  stress  ^  equals  the  sum  of  the  aerodynamical ly 
applied  pressure  -fA  and  any  additional  non- aerodynamic 
loading  -ft?  .  This  then  introduces  the  additional  non-dimen¬ 
sional  parameter. 


<rh 


ft 

<T0 


<rs 


(4.6) 


which  is  to  be  added  to  those  originally  given  in  Eq.  (4.1). 
Summarizing  then  for  a  material  whose  plastic  stress-strain 
curve  is  matched  by  the  Ramberg- Osgood  approximation  as 
stated  previously,  the  general  parameters  for  stress  and 
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deflection  similarity  are, 


(To  j  Vo  )  (Ty  j  c/>q~7o  j  'T^  j 
Eo  0r7)o 

(4.7) 

PbLZ  f  )  ;  E  ,  *,  ^ 

tfo  i-o  Vo  Vo 

The  form  of  the  boundary  condition  Eq.  (4.6)  as 
given  above ,  permits  the  addition  of  non- aerodynamic  loads  ■^r 
on  a  model  to  adjust  for  a  purely  aerodynamic  loading  on 
the  prototype,  provided  this  aerodynamic  loading  is  known 
(see  discussion  in  Section  5.1e,  also  see  Ref.  1).  If  the 
aerodynamic  loading  is  not  known  in  advance,  as  in  aero- 
elastic  type  models,  it  is  then  required  to  simulate  both 
/  <T0  and  separately. 

For  example,  for  such  aerodynamic  type  testing 
over  the  surface  of  bodies  in  hypersonic  flow,  and 

-fn/To  becomes  (see  Eq.  2.3), 


?  (4-8) 

Vo  (7 0 

Maintaining  the  required  aerodynamic  similarity 
for  this  general  body  as  discussed  in  Section  2.2,  the 
above  reduces  to  the  condition  of 

^  (4.9) 

% 
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i  a 


which  replaces  the  (f*  condition  of  Eq.  (4.7)  for 

these  purely  aerodynamic  type  loadings.  The  similarity 
parameters  for  the  combined  aerodynamic,  heat  transfer,  and 
stress  and  deflection  similarity  will  be  discussed  further 
in  Section  5. 

b)  Stresses  and  Deflections  of  Plates 

A  specific  type  of  structure  that  is  often  en¬ 
countered  in  aeroelastic  work  is  the  plate.  Low  aspect 
ratio  wings  are  characteristically  thin  and  deform  much  like 
a  plate.  A  convenient  means  of  representing  the  deflections 
of  a  wing  of  this  type  is  to  transform  the  actual  structure 
into  an  equivalent  plate  of  varying  bending  and  extensional 
stiffness.  It  might  be  of  interest  therefore  to  examine 
the  specific  similarity  parameters  associated  with  heated 
and  laterally  loaded  plates.  See  also  Section  4.3  of  Ref.  1 
where  this  subject  was  also  discussed  along  somewhat  similar 
lines  as  given  below. 

The  governing  differential  equations  of  a 
heated  isotropic  plate  are. 
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where  4  T^y>  ¥>  ~^)  is  the  distributed  lateral  loading  per 
unit  area  of  the  plate. 

These  two  equations  are  the  Von  Karman  equations 
for  large  deflections  of  a  plate  with  variable  stiffness 
properties  under  applied  loads  and  thermal  stresses.  The 
quantity  w(x,y,t)  is  the  lateral  deflection  of  the  neutral 
surface  from  its  initial  position,  F  is  the  stress  function, 
C  and  D  are  the  extensional  and  bending  stiffness  per  unit 
length  respectively,  and  N™  and  Hj,  are  the  thermal  force  and 
moment  per  unit  length  respectively.  These  latter  are 
defined  more  generally  as. 


f/* 

E  J *- 

c  =  /  7^-  > 

J-S/i 
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f  E  <*AT 

Nt=  i—y 
J-Hu 


(4.12) 
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where  °  is  the  plate  thickness  and  AT  is  the  temper¬ 
ature  rise  from  some  initial  value. 

Introducing  the  following  non-dimensional 
parameters. 
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(4.13) 


(4.14) 


The  first  non-dimensional  parameter  appearing  above, 

L*A-f  /p0  Uo  may  be  split  up  further  by  noting  that  the 
lateral  loading  may  be  due  to  aerodynamic,  inertial, 
gravity,  or  externally  applied  forces 
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(4.16) 
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Hence  the  resulting  non-dimensional  parameters  for  a  plate 
under  lateral  loads  and  varying  temperature  are 


L* 'Fa  )  /rrfoL^  n 'rtofrL- 
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(4.17) 


For  the  study  of  small  deflections,  w,  of  heated 
plates,  the  nonlinear  terms  u//d'ix)  appearing 

in  Eq.  (4.11)  may  be  neglected.  In  this  case  the  Uo/l. 
requirement  in  Eq.  (4.15)  would  not  appear.  Further,  it  is 
seen  that  this  would  then  uncouple  the  two  equations  and 
permit  one  to  solve  the  second  equation  for  the  stress 
function  in  the  form. 


F  -  *2°.  f  (  x,  Q)  (4.18) 

Co 

where  f  (*■)?)  is  some  non-dimensional  function  depending  on 
the  temperature  distribution  and  boundary  conditions  of  the 
plate.  Placing  this  back  into  the  first  equation  for  w 
would  result  in  the  following  non-dimensional  parameters  for 
the  small  deflections  of  a  heated,  laterally  loaded  plate 
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4.2  Long-Time  Loading 

In  this  category  the  duration  of  loading  is 
sufficiently  long  that  time  under  load  will  appear  as  an 
additional  parameter  relating  to  such  items  as  creep  and 
fatigue.  The  consideration  of  time  under  load  as  a 
parameter  is  in  addition  to  requiring  the  reproduction  of 
the  parameters  of  Eq.  (4.1)  and  Eq.  (4.6)  for  similitude. 
Reproduction  of  the  actual  time  experienced  by  a  prototype 
for  a  given  loading  will  in  most  cases  be  impossible  be¬ 
cause  of  the  short  operating  times  of  wind  tunnels.  Further¬ 
more,  if  the  actual  time  spent  at  temperature  is  important 
for  a  transient  temperature  test,  this  too  is  impractical  to 
reproduce  since  simulation  of  the  transient  heating  problem 
(for  same  materials  for  example,  Eq.  3.6)  requires  the 
time  scale  to  be  reduced  according  to 


(to)" 


(4.20) 


For  a  steady-state  heating  experiment  this  last  restriction 
however  does  not  apply. 

Reproduction  of  strain-rate  is  important  since  it 
affects  the  creep  and  fatigue  behavior  of  a  structure. 
Therefore  it  is  important  to  consider  both  the  load  level  and 
the  rate  of  loading  of  a  heated  structure  along  with  the 
thermal  loading. 

Duplication  of  the  aforementioned  long-time 
loading  effects  on  a  one-to-one  time  basis  is  impossible 
in  practically  all  cases  if  wind  tunnel  testing  is  utilized. 
Certain  alternatives  are  available  however.  It  is  possible. 
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at  least  for  structures  in  which  a  two  dimensional  state 
of  stress  exists,  to  utilize  results  of  phenomena  such  as 
creep,  tension  and  compression  loadings  at  short  periods 
of  time  and  given  stress  and  temperature  levels  to  extra¬ 
polate  to  longer  periods  of  time  and  different  stress  and 
temperature  levels  with  the  Larson-Miller  parameter 

(4.21) 

where  C  is  a  constant  characteristic  of  the  material  and  t 
is  the  time  to  rupture  or  to  achieve  a  given  strain  (Ref.  27). 
.This  parameter  has  correlated  data  over  a  wide  range  of  times 
and  temperatures  for  many  steels,  high  temperature  alloys, 
and  aluminium  alloys. 

Another  alternative  is  to  utilize  short-time 
experimental  data  in  a  cummulative  sense.  As  suggested  by 
Gerard  (Ref.  28)  for  creep,  the  total  cummulative  creep  may 
be  taken  as  the  sum  of  the  creep  increments  due  to  the  time 
spent  at  each  stress  and  temperature  level.  The  creep  defor¬ 
mation  of  a  prototype,  to  a  first  approximation,  may  thus 
be  accounted  for  by  conducting  several  model  tests  at  a 
given  load  and  temperature  level  where  the  total  testing 
time  is  equal  to  the  prototype  time. 

The  problem  of  fatigue  at  present  is  not  well 
understood,  especially  at  elevated  temperatures.  However, 
some  current  experimental  techniques  can  be  suggested  here. 

The  cumulative-damage  theory  of  Miner  (Ref.  29)  may  be  used 
to  correlate  fatigue  data  for  several  tests  each  run  for  a 
different  number  of  cycles  at  a  specified  stress  and 

temperature  level  with  which  a  fatigue  life  N  is  associated. 
According  to  this  theory,  the  fatigue  life  of  a  prototype 
will  be  reproduced  experimentally  if 
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where  each  term  on  the  left  side  above  represents  a  model 
test  and  each  term  on  the  right  represents  the  percentage 
of  time  spent  by  the  prototype  at  any  given  stress  level. 

The  fatigue  life  N  at  each  stress  level  would  have  to  be 
obtained  from  an  <T-  /V  diagram  at  the  corresponding  temper¬ 
ature  level.  The  cumulative  damage  concept  is  based  on  the 
assumption  that  the  cycling  tests  may  be  conducted  in  any 
desireable  sequence  with  no  effect  on  the  fatigue  life. 

An  alternative  method  consists  of  determining  an 
equivalent  stress  <Tr  of  constant  amplitude  which  would 
have  the  same  overall  effect  as  a  variable  loading  at  different 
stress  levels  <r,  ,  <T,  ,  -  -  ,  and  at  different  cycles  n,  ,  I,.--, 

(Ref.  29).  This  equivalent  stress  (Tv  may  be  expressed  as 


_  (T,*n,  -b  ffl  +  -  -  ~  (4.23) 

n,  +  Oi  -b  -  —  - 

where  the  exponent  x  is  the  inverse  slope  of  the  <T-  N 
diagram  on  a  log-log  plot.  The  equivalent  stress  ^  of 
the  prototype  may  be  estimated  and  then  a  series  of  tests 
may  be  conducted  on  the  model  at  various  stress  levels  and 
cycles  satisfying  Eq.  (4.23)  and  the  condition  (Vv-)„  -  (*)„. 

Additional  considerations  are  important  when 
dealing  with  fatigue.  Scale  effect  may  be  very  significant 
but  the  extent  of  this  has  not  been  fully  understood..  One 
theory  holds  that  a  specimen  in  which  a  large  volume  is 
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subjected  to  high  stress  has  more  chance  of  containing 
a  weak  spot  that  determines  failure  than  does  a  smaller 
piece  (Ref.  30) .  Reproduction  of  geometric  details  and 
use  of  the  same  material  is  also  necessary. 

For  the  problems  discussed  in  this  category  it 
appears  desirable  that  wherever  possible  the  applied  load 
and  temperature  spectrums  of  a  model  should  be  the  same  as 
the  prototype.  In  many  cases  this  will  not  be  possible 
and  some  of  the  suggestions  mentioned  above  may  be  applied, 
accounting  for  these  phenomena  at  least  in  an  approximate 
manner . 
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5.  SIMILARITY  REQUIREMENTS  FOR  COMBINED  AERODYNAMIC, 

HEAT  CONDUCTION,  AND  STRESS -DEFLECTION  STUDIES 

In  the  previous  sections  the  similarity  laws 
required  for  studying  in  detail  the  aerodynamic  pressure 
distribution,  aerodynamic  heat  transfer,  structural  heat 
transfer,  and  stress-deflection  phenomena  at  hypersonic  Mach 
numbers  have  been  discussed  individually.  It  was  indicated 
there  what  similarity  parameters  can  be  expected  in  a 
particular  "restricted  purpose"  test.  Thus,  for  example, 
the  parameters  of  Eq.  (2.10)  are  sufficient  for  studying 
the  inviscid  flow  field  about  rigid  sharp-nosed  slender 
bodies,  the  parameters  of  Eq.  (2.14)  for  viscous  interaction 
effects,  either  Eq.  (2.31)  or  (2.33)  for  studying  the  heat 
transfer  rate  at  the  stagnation  point,  Eq.  (3.40)  for 
ascertaining  the  heat  blocking  effectiveness  of  ablative 
materials  and  Eq.  (4.1)  for  short-time  stress-deflection 
behavior . 

In  the  present  section,  the  combined  similarity  of 
these  phenomena  as  desired  for  complete  aerothermoelastic 
testing  will  be  discussed.  Detailed  combined  simulation 
will  not  be  possible  in  most  cases.  It  is  hoped  however,  to 
bring  out  the  more  significant  parameters  for  these  tests, 
and  to  illustrate  their  application  to  boost  glide  and 
ballistic  type  vehicles. 

5.1  General  Features 

The  similarity  parameters  required  generally  for 
a  complete  aerothermoelastic  model  placed  in  a  perfect  gas 
hypersonic  flow  may  be  summarized  from  Eqs.  (2.5),  (2.28), 
(3.32),  (4.1)  and  (4.6)  as 
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The  model  should  possess  complete  geometric  similarity. 

The  significance  of  each  of  the  above  parameters  has 
already  been  discussed  in  Sections  2,3  and  4.  If  the  above 
similarity  conditions  are  all  met,  the  following  non-dimen¬ 
sional  quantities  will  be  the  same  on  both  model  and  prototype. 


A  u2-  ° 


V~  t  ^(VELoei-rr)  f  U-  (p/iriAC&merrJ  (5.2) 

f~o  to  U 


and  measurements  on  the  model  can  infer  behavior  of  the 
full  scale  prototype.  The  arbitrary  reference  values 
To  Co  to  and  Uo  may  differ  on  both  model  and  prototype. 
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but  must  be  consistent  with  the  4th,  5th,  11th,  29th,  30th, 
and  the  32nd  similarity  conditions  of  Eq.  (5.1),  respectively. 

The  parameters  of  Eq.  (5.1)  are  to  be  compared 
with  those  of  the  lower  supersonic  range  given  in  Section  2.4 
of  Ref.  1.  The  absence  of  the  Mach  number  condition  in 
Eq.  (5.1)  is  to  be  noted  for  these  hypersonic  speeds. 

The  parameters  of  Eq.  (5.1)  apply  aerodynamically 
to  a  perfect  gas.  As  mentioned  in  Section  2.1,  because  of 
the  high  temperatures  involved,  real  gas  effects  come  into 
play  and  in  order  to  simulate  these,  the  same  temperatures 
and  same  pressures  in  the  same  gas  should  be  maintained  on 
model  and  prototype.  This  then  requires  adding  the  para¬ 
meters 


U}  }  same  gas  (5.3) 

to  the  list  of  Eq.  (5.1).  Further,  to  account  for  the 
variation  in  material  properties  over  these  wide  temperature 

»  jT"  - 

ranges  (such  items  as  Kg  t  Cf  ,  E  shape  of  stress-strain 
curve,  etc.)  it  would  be  desirable  also  to  use  the  same 
material  for  both  model  and  prototype.  For  similar  temper¬ 
atures,  pressures  and  materials  then,  the  parameters  of 
Eqs.  (5.3)  and  (5.1)  reduce  to 


^  /U  ,  same  §as  ,  Res  =L  ,  to.  ,  €aL 

(5.4) 

'  L>  L"-,  VoL  ,  u;_  ,  fL,  ra  ,  fF 

A  fundamental  conflict  appears  here  in  the  Reynolds  number 
Re  ,  condition.  This  cannot  be  satisfied  now  for  anything 
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other  than  unity  scale  ratio  .  It  appears  then  that 
general  similarity  cannot  be  achieved,  and  that  recourse 
must  be  made  to  less  restrictive  specialized  situations. 

a)  Blunt-Nosed  Body 

Consider  first  the  aerothermoelastic  similarity 
in  the  neighborhood  of  the  nose  of  a  blunt  body.  From 
Eqs.  (2.6),  (2.31),  (3.32),  (4.1)  and  (4.6),  the  similarity 

requirements  are  effectively 
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Here,  the  first  combined  parameter  above  replaces  eight 
of  the  first  nine  parameters  of  Eq.  (5.1).  The  {  ]  and 

["  ]  quantities  do  not  vary  much,  and  it  is  assumed  that  Twt<.7s  • 

*The  other  L  conditions  in  the  8th  and  9th  parameters  are 
not  as  important,  as  they  stem  from  the  detailed  ablation 
flow  requirements.  Also  the  £°/t-  and  ^-o/lO  conditions  can 
be  somewhat  circumvented  by  dealing  with  two  different 
time  scales,  one  for  unsteady  aerodynamic  phenomena  and 
another  for  heat  conduction  phenomena.  These  are  usually 
fairly  well  separated  and  do  not  cause  much  coupling. 
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There  is  now  no  aerodynamic  requirement  of  maintaining 
same  temperatures  and  same  pressures  in  the  gas.  However, 
to  account  for  variation  in  material  properties  and  for 
ablation,  it  would  be  desirable  to  utilize  the  same  materials 
at  the  same  temperatures.  Under  these  conditions,  the 
requirements  of  Eq.  (5.5)  reduce  to 
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now  be  made  compatible  with  the  ninth  parameter  (  U 
requirement),  thereby  permitting  similarity  for  a  non-ablating 
body.  Note  that  if  the  stagnation  enthalpy  is  due  to  velocity 
alone,  ~  *•  Two  of  the  ablation  parameters  (the 

fifth  and  sixth  above)  cannot  be  met  however,  thereby 
excluding  detailed  ablation  similarity  for  this  case.  As 
can  be  seen  from  Section  3.2,  however,  these  two  ablation 
parameters  relate  to  the  flow  of  the  melted  layer.  If  the 
ablation  process  is  such  that  the  melted  layer  is  removed 
very  quickly,  these  parameters  will  not  play  much  of  a  role 
and  hence  may  be  neglected.  Also,  the  unsteady  aerodynamic 
parameter  Ui0/L  might  be  difficult  to  reconcile  with  the 
Io/l  condition  associated  with  the  structural  vibrations 
in  a  dynamic  aeroelastic  test. 

In  shock  tube,  the  stagnation  enthalpy  ^  /Z- 
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b)  Non-ablating  Slender  Bodies  (same  r  ) 

For  the  aerothermoelastic  similarity  of  slender 
bodies  of  the  same  thickness  ratio  'Y  >  the  similarity 


parameters 

as  obtained 

from 

Eqs 

(4.6)  and 

(4.9) 
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The  first  parameter  can  be  dropped  if  for  both 

model  and  prototype.  There  is  now  only  a  weak  aerodynamic 
requirement  to  maintain  same  temperatures  and  pressures  m 
the  gas  stemming  from  viscous  boundary  layer  interaction 
effects  (the  fifth,  sixth,  and  seventh  parameters  above) . 
However,  to  account  for  variation  in  material  properties, 
it  would  be  desirable  to  keep  the  temperatures  the  same 
on  model  and  prototype.  Under  these  conditions,  the 


12),  (4.1) 
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parameters  of  Eq.  (5.7)  reduce  to 


MjAi.  ,  Jt.ULfi,  Jt  ,  e.L  t  A. 

fc  ^  ~  >  f  1  '  f°  '  'fF 


(5.8) 


The  above  are  somewhat  similar  to  the  requirements  for  the 
blunt  body,  Eq.  (5.6).  The  7s  condition  above,  if  adhered 
to  would  generally  require  similar  velocities  U  .  For 
Tw  4<  Ts  and  for  small  viscous  interaction  effects,  this 
Ts  condition  along  with  the  Jfa, UL Jfo  and  Py  conditions 
can  be  disregarded,  thereby  allowing  more  freedom  to  satisfy 
the  parameters  of  Eq.  (5.8)  much  in  the  same  way  as  those  of 
the  blunt  body,  Eq.  (5.6). 

The  parameters  of  Eq.  (5.7),  although  obtained 
specifically  for  slender  bodies,  apply  also  for  arbitrary 
bodies  if  the  Prandtl  boundary  layer  assumption  is  reason¬ 
able  and  there  is  not  much  viscous  interaction.  For  such 
bodies,  the  viscous  interaction  parameters  (fourth  to  sixth 
above)  can  be  disregarded,  and  the  pressure  distribution  is 
controlled  essentially  by  A7,*,  ,  u»/l  .  The  aerodynamic 

heating  is  essentially  controlled  by  the  {ks (^*/tfsTo) 
condition.  See  for  example.  Section  2.3b.  The  heating  here 
has  been  assumed  laminar.  If  it  is  turbulent,  Eq.  (2.35b) 
would  be  used  in  place  of  Eq.  (2.35a)  above.  Again,  the 
unsteady  aerodynamic  parameter  UIo/l  might  be  difficult  to 
reconcile  with  the  condition  associated  with  the 

structural  vibrations  in  a  dynamic  aeroelastic  test. 
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c)  Non-ablating  Slender  Bodies  (different  'Y  ) 

For  a  family  of  affinely  related  slender  lifting 
surfaces  of  differing  thickness  ratios,  aerothermoelastic 
similarity  can  be  achieved  under  certain  conditions.  It 
will  be  assumed  that, 

1)  Model  and  prototype  are  affinely  related  thin 
solid  plates  of  differing  thickness  ratios  'Y  . 

2)  A  flutter  or  static  aeroelastic  type  test  is 
envisaged . 

3)  Small  deflection  plate  theory  is  applicable. 

4)  Heat  flow  into  the  thin  plate  causes  a  very 
small  temperature  variation  through  the  thickness 
of  the  plate. 

Under  these  conditions,  the  similarity  parameters 
are  obtained  from  Eqs.  (2.14),  (2.18),  (2.36),  (3.12),  (4.19) 
and  (2.8),  as 
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For  a  thin  solid  plate 
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where  J  is  simulated  by  the  condition  of  affinely  related 
bodies.  Placing  the  above  into  the  parameters  of  Eq.  (5.9) 
and  rearranging  gives 
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The  A 1*7  condition  may  be  dropped  if  \  .  The 

fifth,  sixth,  and  seventh  parameters  above  refer  to  viscous 
boundary  layer  interaction  effects  and  may  be  dropped  if 
these  effects  are  small. 
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is  desired  to  keep  the  same  temperatures 
in  both  model  and  prototype,  tb~  thermal  stress  condition 
rt'oTo/r1'  would  require  the  same  thickness  ratio  'Y  for 
model  and  prototype  made  of  the  same  material.  The  para¬ 
meters  of  Eq.  (5.11)  are  then  seen  to  reduce  to  those  of 
the  general  body  given  previously  by  Eq.  (5.7). 

d)  Free  Molecule  Flow 

For  a  blunt  body  in  high  speed  free  molecule  flow, 
the  aerothermoelastic  similarity  parameters  are  obtained  from 
Eqs .  (2.23),  (2.45),  (2.49),  (3.6),  (4.1),  and  (4.6),  namely. 
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In  the  above  case,  it  is  not  likely  that  very  high  tempera¬ 
tures  or  plastic  stresses  will  be  reached.  Hence,  working  with 
the  same  material  at  the  same  temperatures  is  not  necessary 
here,  and  some  flexibility  may  be  achieved  by  using  different 
materials  and  different  reference  temperatures  T,  . 

For  slender  bodies,  of  the  same  thickness  ratio  > 

the  corresponding  parameters  as  obtained  from  Eqs.  (2.24), 

(2.47),  (2.49),  (3.6),  (4.1)  and  (4.6)  are  the  same  as  Eq.  (5.12) 

above  plus  the  additional  parameters 


(5.13) 
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For  very  high  speeds  such  that  St  >>  X  ,  the  above  para¬ 
meters  drop  out  leaving  only  those  of  Eq.  (5.12)  again. 

For  affinely  related  slender  bodies  of  differing 
thickness  ratio  r  ,  aerothermoelastic  similarity  can  be 
achieved  subject  to  the  assumptions  mentioned  in  Section  5.1c. 
Under  these  conditions,  the  similarity  parameters  as  obtained 
from  Eqs.  (2.24),  (2.47),  (2.49),  (3.12),  (4.19),  and  (2.21) 
are 
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the  first  three  parameters  would  drop 

out.  Also,  for 

the  same  temperatures  in  both  model 

and 

prototype  material,  the  J.To/r2-  condition  requires  the 
same  'Y  ,  and  the  above  reduces  to  the  conditions  of  Eq.  (5.12). 
It  is  to  be  noted,  in  free  molecule  flow,  the  temperatures 
are  not  apt  to  become  very  high.  It  is  then  not  so  necessary 
to  reproduce  the  same  materials  and  the  same  temperatures 
as  it  would  be  if  higher  operating  temper*atures  are  present. 

e)  Incomplete  Aerothermoelastic  Testing 

By  incomplete  aerothermoelastic  testing  is  meant 
here  the  use  of  other  auxiliary  means  than  just  the  air 
stream  alone  to  provide  the  net  aerodynamic  pressure  input 
-fver  and  the  net  heat  transfer  input  ^g/er  to  the  model 
structure.  This  requires  a  priori  knowledge  of  what  these 
loadings  -fue-r  and  will  be  on  the  prototype.  In  all 


70 


the  preceding  similarity  investigations  of  Section  5,  a 
complete  aerothermoelastic  model  was  envisaged,  i.e.,  one 
which  when  immersed  in  a  high  speed  flow  would  simulate 
automatically  the  proper  aerodynamic  pressure  and  aerodynamic 
heating  environment  of  the  prototype.  Since  as  was  seen  it 
is  difficult  to  obtain  complete  aerothermoelastic  models 
except  in  certain  specialized  situations,  recourse  can  be  had 
as  a  last  resort  to  the  cruder  incomplete  aerothermoelastic 
testing. 

Looking  at  the  problem  of  a  body  subjected  to  a 
general  force  distribution  anc^  a  general  heating  rate 

prr  ,  the  similarity  parameters  may  be  summarized  from 


Eqs.  (3.10), 

(3.30), 

(4.1),  and  (4.6)  as 
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In  the  above. 

the  f i: 

rst  four  parameters 

stem  from  the  heat 

conduction  requirements,  the  fifth  to  fifteenth  from  the 
ablation  requirements,  and  the  remaining  from  the  stress 
and  deflection  requirements.  The  net  inputs  to  the  system 
are  defined  as 
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where  and  are  any  additional  pressures  and  heating 

rates  to  be  applied  to  the  model  in  order  to  satisfy  the 
basic  input  parameters  above, 

±  (5.17) 

Eo  KoTo 

on  both  model  and  prototype.  In  evaluating  'tfJer  and 
estimates  must  be  made  of  the  expected  values  of  aerodynamic 
pressure  ,  aerodynamic  heating  rate  y.A  ,  radiation  loss 

€u/  f  7u/  ,  and  heat  blockage  due  to  vapor  injection 

for  both  prototype  and  model.  It  is  for  this  reason  that  the 
above  testing  procedure  is  called  incomplete  aerothermoelastic 
testing.  Despite  these  limitations,  this  procedure  is  much 
used  for  structural  testing  at  high  temperatures. 

For  the  same  materials  at  the  same  temperatures, 
the  above  parameters  reduce  to 
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If  there  is  no  ablation,  or  if  the  melt  layer  is  removed  very 
rapidly,  the  third,  fourth,  and  fifth  parameters  above  which 
relate  to  flow  of  the  melted  layer  may  be  neglected. 

In  aeroelastic  type  model  tests,  it  is  generally 
difficult  to  add  additional  forces  ^  to  th.e  model.  In 
these  cases,  it  is  expedient  to  simulate  the  'fmtr  completely 
by  aerodynamic  forces  .  As  seen  in  Section  2.2  this 

means  generally  that  =  -fA  -  where  the  non-dimen 

sional  pressure  p  would  be  simulated  by  maintaining  the 
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parameters  of  say  Eqs.  (2.6),  (2.10),  (2.11),  (2.14)  (2.15), 
(2.23),  (2.24),  and  (2.25),  which  ever  is  felt  to  be  more 
appropriate  (usually  for  gross  pressure  effects,  the  viscous 
interaction  and  surface  temperature  effects  can  be  neglected). 
The  conditions  of  Eq.  (5.15)  (say  for  the  same  material  at 
the  same  temperatures)  would  become  then 
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Thus,  in  such  an  aerodynamic  test,  the  dynamic  pressure 
is  essentially  matched.  The  matching  of  the  heat  input  fMerL 
might  be  achieved  by  additional  9F  added  to  the  model,  say 
by  radiant  heaters  in  the  tunnel  walls.  The  accurate 
distribution  of  the  heating  rate  over  the  model  might  present 
some  difficulty.  It  is  to  be  noted  however,  that  adding  heat 
would  probably  be  easier  than  adding  forces  parti¬ 
cularly  in  a  dynamic  aeroelastic  test  where  fF  would  have  to 
vary  largely  and  rapidly  with  time,  whereas  fF  remains  fairly 
constant  with  time. 

5.2  Application  to  Boost  Glide  Vehicles 

The  similarity  requirements  for  two  specific  flight 
trajectories  will  now  be  considered,  namely  the  boost  glide 
vehicle  and  the  ballistic  missile.  The  areas  of  interest 
for  these  vehicles  are  shown  as  regions  II  and  III  on  the 
corridor  plot  of  Fig.  1.  Typical  aerodynamic  conditions  for 
each  vehicle  are  given  in  Fig.  2  (taken  from  Ref.  4),  thereby 
illustrating  the  orders  of  magnitude  of  the  stagnation  temper- 

Other  items  such  as  may  enter  implicitly  through  p  if 

the  Mach  number  is  not  very  high. 

** 

Tests  of  this  nature  have  been  reported  in  Ref.  37,  for  example. 
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ature,  Reynolds  number  per  foot,  Mach  number  and  velocity 
that  may  be  expected  for  each  prototype.  The  significant 
differences between  the  two  trajectories  are  the  higher 
altitude,  and  hence  lower  Reynolds  number,  lower  attendant 
heating  rates,  and  longer  duration  of  heating  for  the  boost 
glide  vehicle  compared  with  the  ballistic  missile. 

The  boost  glide  vehicle,  because  of  its  lower 
aerodynamic  heating  rates  and  longer  duration  of  flight 
time,  would  probably  use  radiation  as  a  principal  means  of 
heat  protection.  It  is  also  apt  to  be  a  vehicle  with  slender 
lifting  surfaces.  The  parameters  of  Eqs.  (5.7)  and  (5.8) 
would  then  be  appropriate.  (It  is  assumed  that  both  model 
and  prototype  operate  at  r  >7  X  . )  Considering  first 

the  case  of  the  same  materials  at  the  same  temperatures, 

Eq.  (5.8)  yields  the  following  primary  significant  parameters 
for  matching. 


r  '  ^  1  ,  €.l  ,  ^ 
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(5.20) 


the  others  being  somewhat  secondary  (viscous  interaction 
effects,  gravity  effects,  Ts  »  7^  ,  etc.).  The  above 
parameters  could  be  matched  except  for  the  radiation 
condition  e0L  ,  and  unsteady  aerodynamic  parameter  Uto/L  > 
by  picking 

'/3 

(fs)r  Lp 
(/<{)*  L-*\  _ 

_  (±)z  <5-21) 


74 


In  the  above,  the  stagnation  enthalpy  ,^s  was  assumed 
caused  by  velocity  alone,  i.e.,  Jts  U*/z.  .  The  two 

characteristic  times  io  from  the  fifth  and  seventh 
conditions  above  are  associated  with  two  different  phenomena, 
heating  and  aeroelasticity  respectively,  and  would  not  gen¬ 
erally  cause  any  coupling  (see  footnote  Section  5.1).  However, 
the  satisfaction  of  the  U*-o/l  condition  from  unsteady  aero¬ 
dynamics  together  with  the  lo/t_  condition  from  the  structural 
dynamics  would  not  be  possible  unless  V  were  the  same.  For 
static  aeroelastic  phenomena  of  course,  both  these  could  be 
neglected.  The  radiation  condition  would  also  cause  some 
difficulty  here  since  it  would  require, 

(€o)a,  =  (5.22) 

(€°)r 

As  the  prototype  usually  operates  at  high  emissivity,  ; 

coating  the  model  with  an  efficient  radiator  (£•)*,  ~  i-  would 
not  achieve  much  in  the  way  of  scale  ratio.  Since  radiation 
is  a  significant  mode  of  heat  transfer  here,  one  must  resort 
to  one  of  the  circumventions  mentioned  in  Section  3.1.  These 
were  to  consider  steady  state  heat  transfer,  to  consider 
different  temperatures  and  materials  for  model  and  prototype, 
or  to  go  to  the  fA/e-r  concept--i.e.,  the  incomplete  aero¬ 
elastic  testing  of  Section  5.1e. 

Utilizing  this  incomplete  aeroelastic  testing 
technique  the  similarity  parameters  at  high  Mach  numbers  for 
the  same  materials  at  the  same  temperatures  become  basically 
(see  Eq.  5.19) 
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where  some  estimate  of  the  net  heating  rate,  ?*er  ~  f/\  ~  f 
+  fF  would  have  to  be  made  for  both  the  prototype 

and  model,  and  the  wind  tunnel  then  would  be  adjusted  to 
conform  with  this.  Additional  heating  by  radiant  heaters 
might  be  used  in  the  tunnel  if  the  tunnel  air  is  not  suffi¬ 
cient  to  provide  the  desired  heating  rates.  The  distribution 
of  the  heating  rate  over  the  model  might  be  difficult  to 
simulate  accurately  by  these  artificial  means.  Notice  that 
this  incomplete  aeroelastic  testing  can  also  satisfy  the 
time  conditions  bo/i_  and  ULo/l  for  dynamic  aeroelastic 
tests  by  having  U  the  same  for  model  and  prototype. 

5.3  Application  to  Ballistic  Missiles 

The  ballistic  missile  is  generally  a  blunt  nosed, 
wingless  vehicle  that  experiences  a  very  high  but  short 
duration  aerodynamic  heating  rate  over  its  trajectory.  This 
heating  occurs  mainly  in  the  vicinity  of  the  stagnation 
point  and  is  apt  to  be  small  elsewhere.  Ablation  is  usually 
the  primary  means  of  thermal  protection  for  this  type  of 
vehicle  and  radiation  does  not  play  a  large  role.  The  para¬ 
meters  of  Eqs .  (5.5)  and  (5.6)  would  then  be  appropriate. 
Considering  first  the  case  of  the  same  materials  at  the 
same  temperatures,  Eq.  (5.6)  yields  the  following  significant 
parameters  for  matching 
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the  others  being  somewhat  secondary  (radiation,  gravity 
effects,  7~s  >>  Tw  ,  etc.).  It  is  immediately  seen  that  there 
is  difficulty  in  reconciling  two  of  the  ablation  parameters 
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(fourth  and  fifth)  with  the  heating  (first)  and  aeroelastic 
(seventh)  requirements  here.  It  appears  then  that  detailed 
ablation  similarity  for  these  cases  cannot  be  maintained. 

As  mentioned  in  Section  3.2,  however,  these  parameters  relate 
to  the  flow  of  the  liquid  layer.  If  this  is  removed  very 
quickly  or  if  the  ablation  is  essentially  a  sublimation  type 
process,  these  parameters  would  not  appear  and  hence  may  be 
neglected.  There  is  also  some  difficulty  in  matching  the 
unsteady  aerodynamic  condition  Ui,/L  together  with  the 
structural  dynamics  condition  without  going  to  the  same 

velocity  U  .  For  static  aeroelastic  type  phenomena,  both 
of  these  conditions  could  be  neglected.  The  similarity  para¬ 
meters  then  appear  to  be  the  same  as  in  the  Boost  Glide  case, 
Eq.  (5.20),  without  the  radiation  term  e,  L  .  The  matching 
of  these  conditions  is  achieved  in  the  same  manner  and  results 
in  Eq.  (5.21)  as  before. 

Alternatively,  one  could  go  to  the  incomplete 
aerothermoelastic  testing  procedure  of  Section  5 . le  where 
r_  •  Application  of  this  results  in  the 

following  basic  similarity  parameters  at  high  Mach  numbers 
for  the  same  materials  at  the  same  temperatures  given  by 
Eq.  (5.19). 
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where  again  some  estimate  of  the  net  heating  rate. 
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would  have  to  be  made,  and  the  tunnel  adjusted  to  conform  to 
this.  The  third,  fourth,  and  fifth  parameters  again  refer 
to  the  melt  layer  and  might  be  neglected  under  appropriate 
circumstances.  Also  the  time  conditions  Uie/L  and  i'/t- 
for  dynamic  aeroelastic  tests  can  be  satisfied  by  having  U 
the  same  for  model  and  prototype. 
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6. 


C ONS IDERATION  OF  TESTING  FACILITIES 


The  parameters  of  the  preceding  sections  indicate 
that  a  hypersonic  wind  tunnel  testing  facility  should  in 
general  have  the  following  capabilities: 

a)  high  free  stream  enthalpy 

b)  high  free  stream  velocity 

c)  high  stagnation  pressure 

d)  large  test  section 

e)  relatively  long  testing  time 

For  a  general  aerothermoelastic  test  duplicating 
the  aerodynamic,  heat  transfer,  and  structural  phenomena 
of  a  prototype  exactly,  it  has  been  shown  that  a  full  scale 
model  must  be  tested  in  a  flow  with  the  same  Ot  and 

-fa,  .  Thus  for  full  scale  similitude,  each  of  the  capabil¬ 
ities  listed  above  is  required.  An  indication  of  the  operating 
wind  tunnel  stagnation  temperatures  and  pressures  (in  the 
reservoir)  necessary  for  duplicating  these  flight  parameters 
(  ?<*>  j  ~fo>  /  ^  i-s  8^ven  Fig.  4  (Ref.  31).  For  example, 

operating  tunnel  stagnation  temperatures  as  high  as  12,000  R 
and  operating  tunnel  stagnation  pressures  as  high  as  10,000  atm 
may  be  required  for  flight  duplication.  It  may  be  noted  from 
this  figure  that  the  tunnel  requirements  for  the  ballistic 
missile,  whose  altitude  range  of  interest  is  below  100,000  ft., 
are  generally  more  severe  than  those  of  the  boost  glide 
vehicle,  which  has  an  altitude  range  of  interest  between 
100,000  and  200,000  ft.  The  high  stagnation  temperatures 
and  pressures  indicated  in  Fig.  4  introduce  practical  engi¬ 
neering  problems  (such  as  melting  of  the  tunnel  nozzle 
throat)  so  that  the  testing  time  at  these  elevated  levels 
must  necessarily  be  short.  Therefore,  wind  tunnel  testing 
times  for  ballistic  trajectory  studies  are  considerably 
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shorter  than  times  for  boost  glide  studies,  placing  a  practical 
limitation  upon  the  characteristic  time  scale  parameter. 

As  an  alternative  to  reproducing  the  flight  parameters 
U >  fao  >  %o  and  in  an  experimental  facility  (which  are  the 

most  stringent  similitude  requirements),  the  possibility  of 
simulating  the  parameters  associated  with  the  specialized 
phenomena  discussed  in  the  preceding  sections  and  of  utilizing 
the  concept  of  incomplete  aerothermoelastic  testing  of  Section 
5.1e  is  of  interest. 

The  pressure  distribution  about  blunt  bodies  in 
hypersonic  flow  (Eq.  2.6)  is  seen  to  depend  primarily  on 
geometric  similarity  (/3  )  provided  the  flow  Mach  number  is 
approximately  3  or  greater.  No  condition  is  imposed  on  Res 
(and  consequently  on  L)  or  on  ambient  free  stream  conditions. 
Therefore  the  pertinent  tunnel  capability  requirement  here 
for  similitude  of  ~T  H*s  is  simply  that  Ma>  be  greater  than 
approximately  3.  Conventional  blow-down  hypersonic  wind 
tunnels  which  can  simulate  up  to  approximately  10  are 

thus  adequate  here.  Gases  other  than  air,  such  as  helium, 
may  also  be  utilized  here  since  the  effect  of  y  is  absorbed 
in  'fs  (Eq.  2.7).  Heat  transfer  studies  about  blunt  bodies 
however  introduce  the  additional  requirement  of  a  high  enthalpy 
airstream  ( J,s  )  along  with  a  consideration  of  Rcs  (Eq.  2.31). 
The  shock  tube,  which  is  capable  of  producing  a  high  enthalpy, 
low  f^<o  flow  has  been  used  extensively  for  this  purpose. 

This  facility  is  capable  of  reproducing  enthalpies  corres¬ 
ponding  to  near  satellite  velocities  (Fig.  5).  Hypersonic 
wind  tunnels  on  the  other  hand  are  limited  in  this  respect 
in  that  they  can  reproduce  flight  enthalpies  corresponding 
to  M*>  «  6  only  with  the  aid  of  special  heaters.  Ordinarily, 
however,  these  tunnels  are  limited  to  stagnation  temperatures 
of  about  1000°F,  corresponding  to  flight  at  approximately 
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/%,  ~3.5 in  the  stratosphere.  The  principal  limitations  of 
the  shock  tube  are  its  extremely  short  testing  time  (under 
1  millisecond)  and  small  test  section,  whereas  the  blowdown 
tunnels  have  large  test  sections  and  times  up  to  about  one 
minute. 

Pressure  distribution  studies  about  slender  bodies 
(Section  2.2)  require  facilities  with  a  higher  Mach  number 
capability  than  is  generally  necessary  for  blunt  body  experi¬ 
ments.  Here  again  hypersonic  wind  tunnels  are  useful  for 
inviscid  pressure  experiments  where  the  pertinent  flow  para¬ 
meters  are  Mr,  ~  ,  y  •  Mach  numbers  up  to  approximately  15 
(using  air)  can  be  reproduced  (Fig.  5 )  at  the  corresponding 
stagnation  temperatures  required  to  prevent  condensation.  Still 
higher  Mach  numbers (/M* a zo )  can  be  achieved  with  helium  tunnels 
where  condensation  is  not  a  problem.  Some  latitude  in  scal¬ 
ing  is  possible  here  since  T  and  J  may  be  adjusted  within 
certain  practical  limits.  Note  that  there  is  no  condition  on 
the  length  scale  L  here  since  there  is  no  requirement  on  Res  . 
Viscous  interaction  effects  and  heat  transfer  studies  intro¬ 
duce  the  additional  considerations  of  stagnation  temperature 
and  Res  (Eqs.  2.14  and  2.37).  The  blow-down  wind  tunnel, 
as  already  indicated,  is  limited  with  respect  to  stagnation 
temperature  but  may  be  used  for  the  lower  velocity  region  of 
the  boost  glide  trajectory,  and  shock  tubes  are  limited  with 
respect  to  Mach  number.  Higher  enthalpies  and  Mach  numbers 
are  possible  in  the  shock  tube  tunnel  (Ref.  31)  which  expands 
the  high  enthalpy  flow  of  the  shock  tube  through  a  nozzle. 

Fig.  6,  for  example,  illustrates  the  operating  range  of  the 
Cornell  Aeronautical  Laboratory's  24  inch  shock  tube  tunnel. 
Mach  numbers  up  to  approximately  20,  stagnation  pressures  up 
to  about  400  atmospheres  and  stagnation  temperatures  as  high 
as  6200°F,  may  be  achieved  in  this  tunnel.  The  shock  tube 

£  ' 

Flight  stagnation  temperatures  will  not  be  reproduced  here, 
however,  for  M  greater  than  about  8. 
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tunnel  thus  has  the  capability  of  reproducing  the  required 
aerodynamic  parameters  associated  with  a  major  portion  of 
both  boost  glide  and  ballistic  missile  trajectories  in  reason¬ 
ably  large  test  sections.  As  with  the  shock  tube,  the  shock 
tube  tunnel  is  limited  with  regard  to  testing  time  (order  of 
5  to  15  milliseconds) . 

Another  facility  which  makes  use  of  shock  wave 
heating  to  produce  high  enthalpy  flow  at  high  and  high 

stagnation  pressures  is  the  Wave  Superheater  being  developed 
at  Cornell  Aeronautical  Laboratory  (Ref.  32).  Mach  numbers 
of  5  to  15  and  air  temperatures  between  5000°R  and  10,000°R 
are  possible  for  flow  durations  up  to  approximately  30  seconds. 

Very  high  stagnation  enthalpies  and  Mach  numbers, 
corresponding  to  ballistic  trajectory  values  are  also  possible 
in  the  arc-discharge  or  "hotshot"  wind  tunnels.  Reynolds 
numbers  are  an  order  of  magnitude  greater  than  those  produced 
by  shock  tunnels  but  running  times  are  only  slightly  greater. 
These  tunnels  are  particularly  suitable  for  high  /V^  ,  high 

enthalpy  flows  and  their  chief  disadvantage,  in  addition  to 
short  running  times,  is  the  contamination  of  the  flow  due 
to  the  arc  burning  process. 

Inevitably,  recourse  must  be  made  to  incomplete 
aerothermoelastic  testing  as  discussed  in  Section  5.1e.  As 
already  indicated  there,  it  is  generally  more  difficult  to 
artificially  apply  the  aerodynamic  pressure  environment 
than  the  heat  transfer  environment,  particular ly  for  dynamic 
aeroelastic  type  tests.  An  estimate  of  each  is  first  neces¬ 
sary.  The  pressure  and  heat  transfer  environments  however 
may  be  adequately  estimated  in  static  type  situations  by 
utilizing  experiments  conducted  in  the  facilities  discussed 
earlier  in  this  section,  such  as  the  shock  tube  tunnel  for 
example,  in  conjunction  with  theoretical  estimates.  With 
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this  information,  tests  which  are  intended  primarily  for  the 
study  of  structural  and  material  behavior  and  which  require 
large  test  sections  and  long  running  times,  may  be  conducted 
in  facilities  such  as  rocket  exhausts,  plasma  jets,  radiant 
heaters,  etc.  These  latter  facilities  are  particularly 
attractive  because  of  the  high  thermal  fluxes  which  they 
can  provide. 

The  development  of  rarefied  gas  flow  facilities 
has  lagged  the  continuum  flow  facilities.  Presently  opera¬ 
ting  low  density  wind  tunnels  are  capable  of  producing 
velocities  of  about  2500  ft/sec.  and  Mach  numbers  up  to 
about  6.  Facilities  are  planned  to  extend  this  range  to 
Mach  numbers  of  approximately  10  (Ref.  33).  Some  typical 
tests  on  cylinders  conducted  in  low  density  wind  tunnels 
are  reported  in  Refs.  13,  11  and  34.  Knudsen  numbers  as 
high  as  185  (Ref.  13)  and  molecular  speed  ratios  of  the  order 
of  2  to  3  were  achieved.  Tests  of  a  different  nature  using 
a  model  supported  by  an  arm  revolving  at  high  speed  in  a 
vacuum  chamber  are  described  in  Ref.  35.  The  pressure  in 
the  chamber  was  varied  and  molecular  speed  ratios  of  approx¬ 
imately  1.5  were  attained. 

The  regions  of  simulation  of  some  typical  aero¬ 
dynamic  facilities  are  indicated  in  Fig.  5  along  with  typical 
boost-glide  and  ballistic  missile  trajectories  (Ref.  33) . 
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7. 


SUMMARY  AND  CONCLUSIONS 


The  present  report  has  derived  the  similarity  para¬ 
meters  required  for  aerothermoelastic  studies  at  Mach  numbers 
extending  into  the  hypersonic  regime.  This  study  has  consid¬ 
ered  in  detail  the  external  pressure  and  heat  flux  imposed 
upon  the  structure,  the  heat  transfer  distribution  within  the 
structure  (for  both  ablating  and  non-ablating  materials)  and 
the  stress  deflection  behavior  under  short  time  and  long  time 
loading  conditions. 

Similitude  for  the  general  hypersonic  aerothermo¬ 
elastic  problem  is  not  possible  and  a  model  must  be  an  identi¬ 
cal  size  replica  of  the  prototype  with  the  same  pressure  and 
thermal  environments.  Although  the  Mach  number  does  not  play 
a  dominant  role  here,  the  basic  conflict  in  similitude  arises 
from  mutually  satisfying  the  requirements  of  same  temperature  T, 
Reynolds  number  f?es  ,  radiation  parameter  £B~n>  JLj/<c  ,  and 
the  aeroelastic  parameter  f^U^/Eo 

Certain  relaxations  of  this  basic  conflict  are 
possible  when  considering  specialized  situations.  This  has 
been  brought  out  in  the  parameters  which  apply  specifically 
to  blunt  bodies  and  affinely  related  slender  bodies.  These 
specialized  situations  were  examined  for  both  continuum  and 
free  molecule  flow.  Even  though  these  specialized  situations 
offer  some  relaxation,  the  problem  of  similitude  is  still 
extremely  difficult  because  of  the  high  operating  tempera¬ 
tures  and  the  attendant  problems  of  creep  and  fatigue. 

Recourse  may  then  be  made  to  "incomplete  aerothermoelastic 
testing"  in  which  the  pressure  and  thermal  loadings  must  be 
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estimated  in  advance  and  applied  artificially  to  the  model. 
Otherwise  it  appears  that  only  ’’restricted  purpose"  models 
investigating  some  specific  problem  such  as  aerodynamic 
pressure,  aerodynamic  heating,  heat  conduction,  or  stress 
and  deflection  behavior,  separately  is  the  only  other 
alternative  here,  rather  than  a  complete  aerothermoelastic 
model . 
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INITIAL  POSITION  BEFORE  ABLATION 
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FIGURE  3.  COORDINATE  SYSTEM  FOR  ABLATING  BODY. 
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FIGURE  6 

OPERATING  RANGE  OF  CORNELL  AERONAUTICAL 
LABORATORY'S  24  INCH  SHOCK  TUNNEL. 

TAKEN  FROM  REF.  36 


97 


